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ABSTRACT 


A MW level FEL is being designed with a short a Rayleigh length resonator to 
increase the spot size at the mirrors and to avoid mirror damage. In this thesis, it is found 
that it is desirable to focus the electron beam to improve the FEL extraction efficiency. 
Three-dimensional simulations show that the focused electron beam increases the 
extraction efficiency far beyond the desired value of 0.7%. It is also found in this thesis 
that shifting the electron beam off-axis less than 0.3 mm, the efficiency remains above 
the required value. The proposed FEL design uses high power, short optical pulses whose 
spectrum may cover many absorption lines. The absorbed laser energy can heat up the air 
resulting in defocusing the laser beam (thermal blooming). This thesis shows that thermal 
blooming is not an issue for a moderate clear atmosphere when the stagnation zone size 
remains less than 10 m. A transitional step for the development of a MW level FEL 
weapon is the proposed 100 kW upgrade of the Thomas Jefferson National Accelerator 
Facility’s FEL. It has also been shown in this thesis that the use of a step-taper undulator 
slightly improves the performance of the FEL. Finally, the potential of various high 


average power solid-state laser designs are reviewed. 
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EXECUTIVE SUMMARY 


At the Naval Postgraduate School, a high power (MW level) FEL for naval 
applications is being designed. The proposed FEL uses a short Rayleigh undulator, which 
increases the spot size of the optical beam at the mirrors and leaves a small optical mode 


waist in the center of the cavity. 


This thesis shows that the small optical mode waist implies that an electron beam 
shift that is induced vibrationally may reduce the overlap between the electrons and the 
optical mode in the interaction region and thus result in lower gain and efficiency. 
Multimode simulations show that the electron beam shift can change the optical 
wavelength. The peak extraction efficiency steadily decreases as the electron beam is 
offset further from the undulator axis. But when the electron beam shift is less than 0.3 
mm, it is found that the extraction efficiency remains above the goal value of 0.7% for 
the 1 MW output power. This value of the electron beam shift is well beyond the 


experimental design tolerance of 0.01 mm shift. 


In this thesis the electron beam focusing was also studied, by varying the electron 
waist radius while keeping the total current fixed. This process improves the overlap 
between the electron beam and the optical mode and increases the extraction efficiency 
far beyond the desired value of 0.7%. This indicates that it may be possible to obtain MW 
operation at a lower average current. The induced electron energy spread is increased 


slightly but remains within the design goal of 15%. 


The proposed FEL uses high power optical micropulses 1 ps long and a pulse 
repetition frequency of 02=750 Mhz. The pulse mode of operation leads to frequency 
spread depending on the FEL pulse shape in the time domain. The spectrum may cover 
many absorption lines. The absorbed laser energy heats the air resulting in local 
modification of the refractive index, and causing defocusing of the laser beam (thermal 
blooming). It is found that thermal blooming is not an issue for a moderately clear 
atmosphere when the stagnation zone size remains less than 10 m. 


XXx1 


A transitional step for the development of a MW level FEL weapon is the 
proposed 100 kW upgrade of the Thomas Jefferson National Accelerator Facility’s FEL. 
It has been shown in this thesis that there is a slight improvement in efficiency with the 
use of a small negative step-tapered undulator, while the induced energy spread stays 


well below the 15% limit. 


The solid-state laser also offers significant potential for a directed energy weapon 
system. This type of laser has already demonstrated a kW level average output power. 
However, in order for solid state lasers to realize their potential for directed energy 
weapons application, technological advances are required in both the supporting 
technologies and the laser itself. In order to make progress, promising laser designs, such 
as the fiber, continuously cooled, and heat capacity lasers, must develop independent 
technologies dealing with the pump process, cost reduction of the diode pump. Laser 
materials, thermal management, and beam control will play critical roles in the weapon 


design. 
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I. INTRODUCTION 


The requirement for supporting new warfare concepts over a wide range of 
missions has prioritized the development and availability of directed-energy weapons 
(DEW). The potential for a speed-of- light response can add new dimensions and improve 
numerous aspects of military operations, from initial target detection and identification, 
to battle damage assessment after the attack. DEW advantages include speed-of-light 
engagement, high-probability hard kill with unique damage mechanisms and multiple 


target engagement capability. 


An example of a DEW is a high-energy laser. High-energy lasers (HELs) are 
extremely fast and extremely precise. They begin their attack within seconds of the 
target’s detection and complete its destruction a few seconds later. This enables them to 
execute multiple shots in order to destroy the target or engage multiple targets. In 
addition, HELs can engage with smaller bursts of energy as needed for nor lethal attacks, 
offering the potential of an intermediate response option as a complement to existing 


weapon systems. 


In principle, these weapons can be used in either offensive or defensive 
applications, though the predominant role has been defensive. The potential missions 
include ballistic missile defense, air defense, attack against ground and maritime targets, 
maritime self-defense, space control, urban operations, clearance of exposed mines, 
surveillance and counter-surveillance. The platforms can be ships, large or tactical 


aircraft including helicopters, ground vehicles, ground bases and spacecraft. 


A laser beam delivers its energy either in a continuous wave mode (CW) or in a 
pulsed mode to a relatively small spot on the target. The incident intensity can be 
sufficient to melt through a missile body in a few seconds. Moreover, if the heated area is 
under stress (aerodynamic or static pressure load), the lethal effects to the target system 
occur more quickly. The beam aims at specific points of the target’s body that are known 


to be vulnerable such as fuel tanks and aerodynamic control surfaces of a missile. 


A HEL weapon system requires more than just a laser device with a specific 
power level. It also requires command, control, communication and intelligence 
subsystems, which control and monitor the weapon. The system should be able to provide 
battle management, including target acquisition, engagement control, kill assessment 


communication and compatibility with other assets or weapon systems of the carrier unit. 


The laser device produces high-energy light at a particular wavelength or small 
set of wavelengths. Different device technologies will produce widely different power 
levels and beam quality, which is a measure of how effective the laser is in putting its 
light into a desired spot size at the target. For HEL weapons, good beam quality is as 
important as high power. Once the laser has produced high power, the energy must be 
transmitted to the target. The beam director directs the laser power toward the target in 
such a way that the laser beam can deliver a lethal intensity on the target. A schematic of 


the system integration into the host platform of a ship is shown in Figure 1. 


Beam Director 
Ship Deck 


Command 
Control 
Communication 





Figure 1. Schematic of a Ship-Based HEL Weapon System. (After [Ref. 1]). 


When a laser beam propagates through the atmosphere, molecular constituents 
and air density fluctuations cause scatter and absorption. Absorption can initiate thermal 
blooming if the power density is high enough, and air density and thermal gradients can 
cause turbulence. These effects may attenuate and spread the beam. The effects are 


typically worse near the surface of the earth than at high altitudes. Figure 2 shows an 


example of a maritime environment in which the high-energy laser weapon system has to 


operate. Nevertheless, three fundamental parameters affect the beam on its way to the 


target: 


Wavelength: Shorter wavelength light is less affected by diffraction but 
can be more strongly affected by device optical imperfections, turbulence 
and other atmospheric effects. 





Figure 2. A Day at Sea. (From [Ref. 1]). 


Transmitter diameter: A large transmitter is able to produce a small spot 
on the target when the beam is diffraction limited, and conversely, a small 
transmitter can produce only a larger spot. Down to a certain target 
dependent limit, small spots are desirable because they have higher power 
density and can be more lethal to a target. 


Range: The distance of the target from the HEL weapon system is a 
critical parameter that strongly influences system design. Destruction of 
distant targets requires much tighter beam control tolerances. 


The three kinds of laser devices that are candidates for HEL weapons are: 


Chemical lasers. Chemical lasers produce radiation from excited atomic 
or molecular species resulting from chemical reactions. Weapon class 
chemical lasers include hydrogen and deuterium fluoride (HF and DF) and 
chemical oxygen/iodine lasers (COIL). These devices have achieved 
megawatt power levels with good beam quality. [Ref. 2, 3]. 


e Solid-state lasers. Solid-state lasers produce radiation from excited 
atomic or ionic species of a solid-state material. The atomic excitation 
succeeds by providing a power supply to the solid-state material. Solid- 
state lasers include heat capacity, fibers and continuously-cooled lasers 
and have produced high average power, at the kilowatt level with good 
beam quality. [Ref. 2, 3]. 


e Free Electron Lasers (FEL). Free electron lasers produce radiation by 
passing a relativistic electron beam through an alternating magnetic field. 
These devices have produced kilowatts of average power with good beam 
quality. [Ref. 2]. 


Currently, each of these device types has limitations, such as the limited magazine size of 
current chemical lasers, the strong thermal management and power supply requirements 


of current solid-state lasers, and the large size and weight of current free electron lasers. 


In the maritime environment, the threats include high speed, high maneuvering 
anti-ship cruise missiles (ASCMs) and fast attack boats. Navy units are currently able to 
engage these threats using kinetic kill systems. The effectiveness of these weapons 
against the potential threat can be severely reduced by the decreased reaction time 
available, especially in littoral operations. A HEL has the potential to engage over a wide 
spectrum of threats with precision, accuracy and almost zero reaction time, as shown in 


Figure 3. 





Figure 3. A Ship in Self- Defense Using HEL Weapon Systems. (From [Ref. 3]) 


Given the electrical power availability on modern navy ships, the FEL and the 
solid state laser are logical candidates for maritime HEL weapon systems, even though 
only the chemical lasers have demonstrated weapornlevel powers. A_ significant 
advantage of FELs and SSLs is that they do not generate the toxic effluents of chemical 
lasers. The issue of handling the chemical laser on board is serious since the crew could 
be harmed and equipment damaged. Another disadvantage is the need for a logistics 
supply of hazardous and toxic chemicals. Furthermore, the chemical lasers are large and 


heavy, which makes them difficult to exist on a compact configuration. 


Chapter II gives an overview of the theoretical background and discusses the 
physics behind the FEL. Starting from the Lorentz force equation the dimensionless 
pendulum equation is obtained, which describes the electron phase space evolution. In 
this way information about the electrons’ energy exchange with the optical wave is 
provided. Additionally, beginning from the Maxwell’s wave equation the optical wave 
equation is obtained. The optical wave equation explains the light wave propagation 


through the undulator. 


Chapter III of this thesis describes the propagation consequences of short pulse 
free electron lasers. Specifically, the propagation issue of thermal blooming related to the 
power spread in the frequency domain that is caused by the pulse format of the FEL 


weapon system will be discussed. 


Chapter IV describes the results of simulations for the proposed 1 MW FEL 
weapon system, using a short Rayleigh length undulator. Three-dimensional simulations 
are used to study the effects of varying the electron beam radius and angular spread to 
enhance FEL gain and efficiency. Additionally, the effects of off-axis shifting of the 


electron beam are studied. 


Chapter V presents the simulation results of the upgraded to 100 kW Thomas 
Jefferson National Accelerator Facility (TJNAF) FEL using a step-tapered undulator. 


Chapter VI discusses the physics behind the solid-state lasers and presents their 


potential to be developed as directed energy weapon systems. 


Chapter VII summarizes the results found in the present study. 


Il. FEL THEORY 


A. INTRODUCTION 

An FEL is a device that extracts kinetic energy from a relativistic electron beam 
and converts it into coherent electromagnetic radiation. The radiation can be over a large 
part of the spectrum, even at wavelengths that conventional sources cannot operate. 
Therefore, the FEL is an important tool for materials science research, chemical 


technology, biophysical science, medical, and military applications. 


Lasing occurs by sending the relativistic free electrons through a periodic 
magnetic field called a wiggler or undulator. The wiggler induces transverse forces on the 
electrons which undergo a wiggling motion as they travel through it. This results in 
electromagnetic radiation emitted in the forward direction of the electron beam. The 
wiggler and the radiation can be combined to produce a beat wave, called a 
ponderomotive wave, that travels slower than the radiation and can be in synchronism 


with the electrons. 


The basic elements of a typical FEL are displayed in Figure 4. The injector 
generates short periodic bursts of electrons from a photocathode, which is illuminated by 
a drive laser. Typical peak injector currents are in the range of | to 100 amps with pulse 
widths ranging from 1 to 50 ps [Ref. 2]. The injected electrons are fed into the 
accelerator, where they are accelerated by an electric field to nearly the speed of light. 


Thus, the electron obtains kinetic energy, 


K,=(y¥- Lm,c’, (2.1) 


1 1 
y= ea (2.2) 


1-(v,/c)’  Vl-B-B 


where y is the Lorentz factor, v 


e 


is the electron velocity, c is the speed of light in 


vacuum, B =Vv,/c is the dimensionless electron velocity, m, is the electron rest mass, 


and m,c” =0.511 MeV is the electron rest energy. Typical accelerators operate in the 


range of 20 MeV to 500 MeV where y ranges from 40 to 1000 [Ref. 2]. The most 
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common electron accelerators are the radio frequency linear accelerator (RF Linac) and 
the superconducting RF Linac. The latter operates at the liquid Helium temperature of 


2°K and is used at the Jefferson Lab FEL. 






Injector 


Figure 4. Major FEL Components. (From [Ref. 4]). 


The next step is to guide the accelerated electrons into the wiggler with the help 
of bending magnets. The interaction between the relativistic electron beam and the 
periodic magnetic field of the wiggler creates and amplifies the optical field. Alternating 
the polarity of a series of permanent magnets or electromagnets produces the periodic 
magnetic field. The wiggler is oriented along the z-axis, with the magnets parallel to the 
y-axis. 

The wiggler is characterized by the dimensionless undulator parameter K given by 
[Ref. 5] 


eB A, 
Oe a 


where A, = L/N is the undulator wavelength, L is the undulator length, N is the number 
of undulator periods, e is the electron charge magnitude, and B,,, is the rms undulator 


field strength. 


Only a small portion of the electron energy is converted to optical energy. When 
the electron beam exits the wiggler it must either be dumped in a radiation shielded dump 
or fed back into the accelerator to return its excess energy back to the system. For the 
latter case, the accelerator is properly phased to decelerate and capture the kinetic energy 
of the electron beam which can be reused to accelerate other electrons. This process is 


called energy recovery and is shown in the schematic of the Jefferson Lab FEL oscillator, 


Figure 5. 


The last major component is the resonator cavity. It is oriented along the wiggler 
axis and consists of an evacuated optical cavity bounded by two mirrors. One of the 
mirrors is partially transmitting to outcouple optical energy from the FEL. The optical 
beam is amplified during successive passes through the wiggler, while a fraction of it is 
allowed to escape. For maximum coupling efficiency between the light and the electrons, 
the optical beam within the resonator should match the electron beam size. 
or 


IR 
wiggler 


inject 
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Figure 5. Jefferson Lab Oscillator FEL Using Energy Recovery. (From [Ref. 4]) 


B. FEL INTERACTION — RESONANCE 


The undulator induces small transverse accelerations of the electrons, which 
results in radiation in the forward direction of the beam. The radiation is a small 


percentage of the electron energy and has the form of a narrow cone of width on the order 
of y'[Ref. 6]. This process is referred to as spontaneous emission and is essential for the 


start-up of an FEL oscillator. 


In the presence of an optical wave, the wiggler mediates the interaction between 


the electrons and the light. The combination of the relativistic electrons, the magnetic 
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field of the wiggler and the optical wave produces a ponderomotive wave. This beat wave 
travels slower than the optical wave and, if it is in resonance with the electrons, there can 
be energy exchange. The ponderomotive wave, which has the same frequency as the 
radiation, causes some electrons to speed up and others to slow down. Thus, the 
ponderomotive wave causes the electrons to bunch. The bunching is critical, as 
stimulated emission occurs when electrons form coherent bunches within an optical 


wavelength. 
At resonance, the optical wave overtakes the electron beam by one optical 
wavelength A while the electron travels over one undulator period 4,. Thus, in a given 


time interval Ar, the light travels a distance 
A+X, =cAt, (2.4) 
while an electron travels 


a, =cB.Ar, (2.5) 


where cB. is the speed that the electron travels along the undulator zaxis. Combining 


the last two equations, the optical wavelength is 


ater, (2.6) 
B. 
Rewriting Equation (2.2), 


1 1 

2 

o> =a a Bes ae (2.7) 
he P= d= $= Gi 

where 8, = (B,. B,.0) is the transverse dimensionless electron velocity component. As 

seen later, the transverse electron velocity can be written 


= 2 


B, =—(ooskyesin by) = Bi= (2.8) 


2 


Combining the last three equations, the optical wavelength becomes 


10 


1+K? 








2 
p=’, (2.9) 

io 

2y 


Since we are dealing with highly relativistic electrons where y >>1 and K is 


usually on the order of unity, the optical wavelength of the FEL at resonance can be 
expressed as a function of the undulator period and the energy of the relativistic electron 


beam as 


2 
od EY (2.10) 
2y Y 








r 


This equation is often referred to as the resonance condition. 
C. ELECTRON DYN AMICS 
1. Pendulum Equation 
This section examines how the ponderomotive wave affects the electron motion 
through the undulator. Relativistic electrons are injected in the z direction and interact 
with a helical undulator magnetic field and the fields of a circular polarized optical plane 
wave, which are described as follows 
B= B(cosk,z,sin k,z,0), 
. = E(siny,cosy,0), (2.11) 
E, = E(cosy,-siny,0), 


1 


where B is the magnetic field amplitude of the undulator, k,=2n/A, is the undulator 


wavenumber, EF is the electric and magnetic field amplitude of the optical wave in cgs 


units, and y is the phase of the optical wave given by 


YW =kz-Or+o, (2.12) 
where k =@/c is the optical wavenumber and 9 is the initial optical phase. 


The Coulomb forces between the electrons can be ignored since the beam energy 
is large. Thus, the total force acting on each electron is approximated by the Lorentz 


force equation [Ref. 6] 
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d(vB Pe oe ee 
ae) =< L 8, +8x(B,+8,)] (2.13) 


Based on the Equation (2.11), the second term in the brackets of the last equation 
can be written 
Bx (B, +B,) = *(-BB. sink,z — EB, cosw )+ 
$ (BB. cosk,z+ EB, cosy )+ (2.14) 
2(BB, sin k,z -—BB,cosk,z+ EB, siny —E B, cosy). 


Then, the x component of Equation (2.13) can be written, 


d(yB.)_ ee 
dt mC 


e 


[ E(1-B. )cosy — BB. sink,z]. (2.15) 


For relativistic electrons, B, ~1 so that E(1—B.)<<BB,. So the last equation 


can be simplified to 


ap. = ee as kz - (2.16) 
dt Yym,c 


Following the same steps for the y component 


dB, eee eBB. 
dt ym,c 





cos kyz. (217) 


The last two equations can be written together as 


i eB aca kz,cosk,z,0), (2.18) 
dt ym,c 
which by integration verifies Equation (2.8), assuming that the electron injection is 


perfect so that the constant of integration is zero. 


From Equation (2.13), with light present, the electron energy exchange is given 
by 


Wy =-*_§.£ =-_£(B, cosy -B, siny), (2.19) 
dt M,C m,C 


and by using Equation (2.18), this becomes 
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we eKE cos(k,z 4y ) = eKE 
m 


m,cy 








cos(C +0), (2.20) 


e 


where C =(k+k,)z—-@r is defined as the electron phase with respect to the 


ponderomotive wave. By differentiating the electron phase twice, a relation for the time 


derivative of the fractional electron velocity in z direction is obtained, 


C =(k+k,)z—-@ =(k+k,) B.c-, (2.21) 
is , fe gH 
C=(k+k,)B.c => babe . (2.22) 


Combining Equations (2.7) and (2.8) and then taking the time derivative on both 


sides, we obtain 





. 2 , 
y> (1+ K?)=1-B? => rate : (2.23) 
Inserting Equation (2.22) into Equation (2.23) we find 
; aie 
Mics BG (2.24) 


Y (1+ K?)(k+k, Je’ 


and since we are dealing with relativistic electrons where B,=1 and k>>k, this 


becomes 


; aie pn 
ej ae a a (2.25) 
Y (1+K*)ke (1+K*)2nc 
This can be simplified further by using Equation (2.10) near resonance so that Equation 
(2.25) becomes 
oa (2.26) 
y 4c 
Finally, combining with Equation (2.22) and solving for the second derivative of the 
electron phase gives 
»  20),eKE 
¢ =—*—cos(€ +0). (2.27) 
m,cY 


e 
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Equation (2.27) has the form of a simple pendulum equation and represents the 


electron phase dynamics. This can be simplified further by introducing the dimensionless 


time t =ct/L and the dimensionless derivatives (...) =d)/dt; 
v=C =|alcos(C +6), (2.28) 


where |a| = 41 NKLE/y?m,c’ is the dimensionless optical field amplitude [Ref. 5] and v 


is the dimensionless electron phase velocity 
v=C =L[(k+k,)B,-k]. (2.29) 


For v=0 we have exact resonance as shown in Equation (2.6). 


At resonance, B. can be written as 


1+ K’ 


=|- 
pals 





(2.30) 


Substituting B. into Equation (2.29) and differentiating, a change in phase velocity Av 


related to a change in electron energy Ay is 
aye ann( (2.31) 
Y 


The FEL pendulum equation is used to describe the microscopic motion of each 


electron over an optical wavelength. The cosine term determines whether each electron 
gains or loses energy. The electrons with phase in the range [—1 /2,n/ 2] absorb energy 
from the optical field and speed up, while those with phase in the range [1/2,3x/2] 


transfer energy to the optical field and slow down. This results in electron bunching in 


each optical wavelength. The other factor of pendulum equation, the dimensionless 


optical field strength, expresses the rate of the electron bunching. For la| >> 7 , the optical 


field is strong and the bunching occurs quickly, while for lal<<a , the optical field is 


weak and the change of the electron phase is small. 
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2. Extraction Efficiency 


The single-pass extraction efficiency n of an FEL can be defined as the fraction 


of the electron beam energy converted to laser light during one passage through the 


undulator, 


y= Sate” __ AY) ba5s 


ym ¥ 
where (Ay ) represents the average energy change of all the electrons after a single pass 
through the undulator. 
From Equations (2.31) and (2.32), the extraction efficiency can be expressed as 
the change of the average electron phase velocity (Av) over a single pass through the 


undulator, 


= _ 4A) . (2.33) 
4n N 


For a typical FEL the change of the average electron phase velocity is (Av) = 2m , so the 


single pass extraction efficiency becomes n = 1/2N . 


D. OPTICAL WAVE EQUATION 
The light wave propagation through the undulator is described by Maxwell’s 


wave equation in the Coulomb gauge, 


Oo Ve . ats 
[v ae =-—J,, (2.34) 


c° Or’ é 


where J, is the transverse current density. The vector potential A is related to the 


circularly polarized plane optical wave of Equation (2.11) by 


B =VxA, (2.35) 
and can be written as a function of z and ¢ as 


A(z,t) = “29 (siny.cosy.0) (2.36) 
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Making the assumption that the amplitude and the phase of the optical field are 


slowly varying both in time and space over an optical wavelength, then 


ae << WE, a << M0, 
Bs : (2.37) 
= << kE, 2 << ko, 


and Equation (2.34) can be simplified to 
(5: 1dE do loo 4nt - 


Set [leo —siny.0)- #( ESF siny.cosy.0) ms ee 


or, 


OE 10E wz 2h FZ -€ 
oz c Ot Ee 





do 106 Dis, oe 
eee ee 
where €,,€, are defined as unit vectors, 
€, =(cosy,—siny ,0), (2.39) 
€, =(siny ,cosy,0). 


The transverse current can be expressed as the sum of each electron charge having 


transverse velocity B,c 


J, =-ecY B.S (x-7) = ec * (cos k,z,sin k, z,0)6° oars (2.40) 
i ey 


where /, is the position of the +th electron. If the number of undulator periods is assumed 


to be large, we can replace the summation over all the electrons with the volume element 


electron density p multiplied by the average phase of the electrons in that volume 


element. Then, from the last two equations, we obtain 


Tf = SAP (cosy +k,z)), 
y (2.41) 
sae 


ol 


6, = 


sin (y + k,z)), 
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where the notation GS corresponds to an average over all the electrons as before. 


Now, combining Equations (2.38) and (2.41) using ¢ +0 =y+k,z gives a 
complex form of the optical wave equation 


d 10 2neKp } _; 
a i a ar) 2.42 
[sts] ‘ y (¢ ) 7 


For the case of a long FEL electron pulse, there is no spatial dependence in the z direction 
and the last equation simplifies to 


1d Be __ 2meKp (e*), 
c Ot y 


which by using dimensionless variables takes the final form of 
> d(lale® 
-_allae®) 


ao j (e* ) (2.43) 


where [a| =4n NKLE/y*m,c’ is the dimensionless optical field amplitude, and 
jJ=8N (en KL) p/y°m,c’ is the dimensionless electron beam current. 


Equation (2.43) is known as the FEL optical wave equation. This equation 
expresses the dependence of the rate of change in the optical field with the dimensionless 
current and the average electron phase changes. 

E. PHASE SPACE ANALYSIS 

The phase space evolution of sampled electrons can be used to describe the 
motion of the electrons within an optical wavelength. A change in electron phase velocity 
is related to a change in electron energy. Therefore, the phase space evolution plot can 
provide information about the electrons’ energy exchange with the optical wave. 

The pendulum equation, which describes the electron phase space evolution, and 
the optical wave equation are coupled together through the dimensionless current density 
j.If j <7, the FEL gain is low, and if 7 >>7 , the FEL gain is high. The phase space 


evolution can be obtained by solving these two equations numerically for a sample 


number of electrons. 
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Suppose a sample of electrons from a monoenergetic electron beam is injected 


into the undulator with a uniformly distributed initial phase ¢, from —2/2 to 3n/2. 
Moreover, the sample electrons are injected at resonance (v, = 0) , and the electron beam 
has low current (j=1) at the upper limit of the weak field regime where a, ~ 2 . Thus, at 
the beginning of the undulator (t =0), the initial conditions for the electron phase C , 


and the phase velocity v, are, 


(2.44) 


where the index i indicates the i-th electron of the total N, sample electrons. Now 
integration of the pendulum Equation (2.28) gives the phase space evolution in (C,v) of 
the sample electrons. 


The electron evolution in the phase-space plane is shown in Figure 6. The 


electrons appear as colored dots and their evolution is denoted by the change of the color 


from yellow (t =0) to red (t =1). For this case, they follow closed orbits inside the 
separatrix, which is defined by v- = 2\al[1 +sin(G, + )| [Ref. 5]. The sample electrons 
with initial phase C, between [1 25%) 2] gain energy from the optical field and move 
up in phase velocity (v>v,), while an equal number of electrons with initial phase 
between [1/2,32/2] transfer some of their energy to the optical field and move down in 


phase velocity (v<v,). This results in a bunching near the relative phase ¢ =2/2 anda 


zero net energy transfer to the optical field. This is also obvious on the top right plot of 


the same figure, which indicates the optical gain G(t) evolution. The gain at time T is 


defined as 


C= (2.45) 

















-1/2 C 3n/2 0 T 1 


Figure 6. Phase-Space Electron Evolution at Resonance. 


where P(t)o|a(t Nh is the optical power at time t, and P(t=0) «<a; is the initial 
optical power. Furthermore, the gain in the weak field, low current regime can be 
expressed as a function of the dimensionless current j and the initial velocity v, [Ref. 5] 


G(t)=j 2—2cos(vt)—v,t sin(v7 ) . 02.46) 


3 
y 


oO 


The last plot of Figure 6 shows the optical phase evolution (t ). 


The final gain spectrum at t =1, during one passage through the undulator for 


j =1, versus the initial velocity v,, is shown in Figure 7. The horizontal axis of v, can 


be considered a function of the electron beam energy Ay , given by Equation (2.31) and 
2 . . 
centered about the resonant energy y = [2 (1 +K | i; 21. | , or a function of the optical 


wavelength AA = AAv,/2nN centered about the resonant wavelength given by Equation 


(2.10). The gain plot indicates that the gain spectrum is anttsymmetric about v, =0 and 
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shows that for a resonant electron beam there is a zero net energy transfer to the optical 


field. 





13% 
0.0 
- 13% 
-12 Vo 12 
Figure 7. Weak Field Low Current Gain Spectrum. 


In addition, the same plot indicates that the maximum gain (13%) for j=1 is achieved 
slightly above resonance at v, = 2.6, while the maximum loss for the optical field occurs 


for v, =—2.6. 


The electron evolution in the phase-space plane for an electron beam injected at 
v, = 2.6 is shown in Figure 8. Now, the electrons above the separatrix follow open orbits. 
At the end of the undulator, the electrons tend to bunch near ¢ =z and the majority of 


them move down in phase velocity losing energy. This results in a net energy transfer to 
the optical field. The top right gain evolution plot shows significant energy exchange to 


the light wave, where the gain has reached the value of 13%. 
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Figure 8. Phase-Space Electron Evolution with Electrons Injected Slightly Above 
Resonance at v, = 2.6. 


F. SHORT PULSES 


The resonance condition states that the optical wave overtakes the electron beam 
by one optical wavelength A over one undulator period 4,. At the end of the undulator, 


the light wave will be ahead of an electron by NA, which is called the slippage distance. 
The pulses are defined as short pulses in the case where the electron pulse length is 


comparable to or less than the slippage distance. 


Due to spontaneous emission, short optical pulses are produced that bounce 
between the mirrors. The FEL interaction takes places between the stored optical pulse 
bouncing between mirrors and every new electron pulse coming into the undulator. The 
optical pulse and the electron pulse may enter the undulator simultaneously at exact 
synchronism. Otherwise, any displacement between these two pulses defines the 


desynchronism d. 


Suppose a short electron pulse, with the pulse length of the order of the slippage 
distance, enters the undulator in exact synchronism with the optical pulse. Since the 
relativistic electrons travel slightly slower than the speed of light, the optical pulse will be 
ahead of the electron pulse by NA at the end of the undulator. Consequently, the light 


amplification occurs mostly for the trailing portion of the optical pulse, which makes the 
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centroid of the optical pulse travel slower than the speed of light. After some oscillations, 
the optical pulse falls behind the electron pulse and no interaction takes place. Finally, the 


optical power starts to decay due to the cavity loses and the FEL becomes unstable. 


This process can be overcome if the optical pulse arrives slightly before the 
electron pulse at the beginning of the undulator. By adjusting the mirrors position using a 
piezoelectric crystal, the path inside the cavity can be reduced by AS. The displacement 
of the mirrors relative to the slippage distance is the desynchronism, defined by 


d =—2AS/Nid . Typical values for peak steady-state power are d ~ 0.01. 


Less power but more stability is obtained by further increasing desynchronism. 


The final optical pulse width will be longer than the electron pulse, which leads to a 
narrower power spectral density. For even larger values of desynchronism (d>>0.1), 


there will not be sufficient overlap between the two pulses and consequently, the FEL 
steady-state power decreases to zero. 
G. DIFFRACTION 

Now suppose that the amplitude of the optical electric field additionally has x and 


y dependence. The optical vector potential of Equation (2.36) can be written as 


E(x,1) 


A(x,1) = 





(siny cosy ,0). (2.47) 


Beginning again with Maxwell’s wave equation in the Coulomb gauge and 
following the same steps as in paragraph D, we obtain 


Ae ik{ 2922 \ mes = -2BeRO (ec) (2.48) 
2 dz c Of y 


es CoO 
where V, Soop te 
ox” oy 
Now, we introduce the dimensionless time, T =2z/ L=ct/ L, and normalize the 


transverse coordinates x and y to the characteristic mode radius, Lh/u . The wave 


Equation (2.48) takes the form of the parabolic wave equation, where the optical wave is 
driven by the dimensionless electron beam current 


Ze 


i=, O .f =i 
[292 +S Ia ye ja—y(e*). (2.49) 


The operator —V7 i /4 describes the diffraction of the optical beam. In the case where the 
optical beam propagates in free space, j =0 and Equation (2.49) becomes the parabolic 
wave equation for free space 

da(x,y,t) ie 

———— = -V' a(x, y,t), 2.50) 

= 4a + ( » ) ( 

where the dimensionless time t =(z/R)=0-—1 is now the time of evolution of the 
wavefront along the range R. Equation (2.50) indicates that the diffraction is important 


when (Ac / 4Ax’) —1,where Ax represents a structure size of the wavefront. Over the 
full range R, At =1, a large size wavefront (Ax >1/2) will not be severely affected by 


diffraction, while a wavefront with a small size (Av < 1/2) will be affected by diffraction. 


1. Rayleigh Length 
A solution of Equation (2.50) determines the fundamental optical mode inside the 


laser cavity. Typically, the fundamental mode has a Gaussian shape and is determined by 


the cavity mirrors. If W, is the radius of the Gaussian optical mode waist, then 


Z,N=1W, 


an (2.51) 
where Z, is the characteristic spread distance called the Rayleigh length. As the beam 


propagates along z, the spot size W(z) becomes larger in accordance with 


W?(z)=W; (2 ] ; (2.52) 


From the last equation, it is obvious that the laser beam spreads out by diffraction 


and the wavefront area of the optical mode is doubled at a distance equal to the Rayleigh 
length. For Z, = L, Equation (2.51) gives the characteristic mode radius JLa/n , which 


is used for normalizing the x and y coordinates. 
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The mode waist is usually located close to the cavity center. According to 
Equation (2.51) by making Z, very small (short Rayleigh length), the mode waist 
becomes small and the beam expands rapidly due to diffraction as it moves to the mirrors. 
However, the small optical beam waist may reduce the overlap between the optical wave 
and the narrow electron beam, which leads to gain reduction. If Z, is too large, however, 
the mode waist becomes large and the beam does not expand as much as before while 
moving to the mirrors. In this case, a large mode waist leads to gain reduction as it 
reduces the overlap between the optical wave and the narrow electron beam in the center 
of the undulator. Thus, the value of the Rayleigh length is critical for gain optimization. 
In order to minimize the optical mode volume around a smaller electron beam, the 
Rayleigh length should be Z, = L| V2 [Ref. 5]. This is the typical FEL case. 

2. Gaussian Beams 

An optical wave with a Gaussian transverse shape everywhere, and also a solution 
of Equation (2.50), is given by [Ref. 7] 


a ( r,t ) — ay | eles oe) 


we) ; (2.53) 


where a, is the amplitude of the field at the beam waist which is located at T=T,,, @ is 


the optical phase that is given by 


(r,t )=—tan™ § a } or Apee) ; (2.54) 


z zo +(t =) 





'O 


and 


w?(t)=1¢4-— 2, (2.55) 





where z,=Z,/L is the normalized Rayleigh length, and r* =x*+ y’ is the transverse 
position related to the normalized x, y coordinates. The beam spreads out by diffraction as 


it moves away from the position of beam waist T,,. From the last equation, it is obvious 


that the area of the beam doubles when Tt =7,, +z, . 
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This Chapter presented an overview of the physics behind the FEL. The equation 
referred to as the resonance condition was first obtained. This equation expresses the 
optical wavelength of the FEL as a function of the undulator period and the energy of the 
relativistic electron beam. The dimensionless pendulum equation, which describes the 
electron phase space evolution and provides information about the electrons’ energy 
exchange with the optical wave, was also obtained. Additionally, the optical wave 


equation was derived. The optical wave equation explains the light wave propagation. 


Based on the physics behind the FEL, the propagation consequences of short 
pulses from free electron lasers will be discussed in the Chapter that follows. 
Specifically, the propagation issue of thermal blooming related to the power spread in the 


frequency domain that is caused by the pulse format of the FEL will be discussed. 
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Il, PROPAGATION CONSEQUENCES OF SHORT PULSES 
FROM FREE ELECTRON LASERS 


A high-energy laser weapon system operating in a maritime environment faces 
several challenges due to atmospheric effects. The propagation of a light wave through 
the atmosphere just above sea level requires a proper selection of wavelength, power 
level, pulse format, and beam control so that a lethal amount of energy can be delivered 


to the target. 


The primary propagation issues that need to be considered are the absorption and 
scattering. At lower power levels, these effects are linear which means that an increase in 
power at the origin results in an increase in power delivered at the target. However, for 
high power levels, absorption leads to nomlinear effects, such as thermal blooming, 
where an increase in power from the laser results in a decreasing intensity at the target. 
Thermal blooming is a result of the atmospheric heating due to absorption, which creates 
a negative lens in the atmosphere and spreads the beam’s core. There is an upper limit on 
power that can be delivered by a laser without reducing the amount of energy on the 
target. This is known as the critical power, Perit. The value of Pcrit depends on the laser- 


operating wavelength. 


The atmospheric transmittance over a wavelength range is a complex problem and 
crucial to directed energy because it determines the choice of a proper wavelength. Figure 
9 demonstrates the dependence of the absorption coefficient on the wavelength. On the 
horizontal and the vertical axis of that figure are illustrated, in logarithmic scale, the 
wavelength and the absorption coefficient respectively. The wavelengths from different 
kinds of lasers are indicated by red lines. The primary windows with low absorption are 
located at 2.2, 1.6, 1.3 microns (um ), and the lowest one is close to 1.0 um (1.029-1.064 


um ). Thus, the thermal blooming induced by absorption is minimal at this wavelength. 


However, the scattering increases with decreasing wavelengths since it is dependent on 
the refractive index and the mean particle size in the atmosphere. The result 6 that the 
total loss due to absorption and scattering is about equal across the four windows. 


However, simulations using the MOLLY code have shown that P,,i: has a value of about 
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10 MW at 1.0 um for distances of up to a few kilometers while at 1.6 pm is in the 
neighborhood of 3 MW [Ref. 1]. Therefore, if the decision is based on operating the laser 
at megawatts level of power in the maritime environment, the 1.0 um wavelength should 
be chosen. Of course, the 1.6 jm wavelength may be preferable since it falls in the 


region of the spectrum considered “eye-safe”’. 


In reality, the atmospheric conditions of a maritime environment are a function of 
a particular area at a particular time of day. Thus, more than one wavelength may be 
proper to overcome the propagation issues due to the variety of the atmospheric 
conditions. This requirement implies that the free electron laser is the most logical 
candidate for maritime high-energy weapons. FELs have the ability to select a 
wavelength of operation in some bandwidth range, which is a significant advantage 
compared to the chemical or solid state lasers operating in the maritime environment. 
Equation (2.10) indicates that the wavelength of the FEL optical beam A is controlled by 


the undulator period 2,, the undulator parameter K, and the Lorentz factor of the 


relativistic electron beam y , 
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Figure 9. Absorption Coefficient Profile for Maritime Environment. (From [Ref. 4]. 
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A. WHY A MW POWER LEVEL OUTPUT IS REQUIRED 


The kill mechanism of a laser is based on delivering a sufficient amount of energy 
to a relatively small spot on the target. The incident intensity of light should be enough to 
melt through a missile body during the engagement time so that the missile will finally 
break up due to aerodynamic stress. Experiments have shown that a flux density of 10 
kW/cn®’ is sufficient to melt typical missile materials in a dwell time of a couple of 
seconds [Ref. 8]. 

To estimate the required laser output power, examine a typical missile 
engagement scenario. In order to destroy an incoming missile, the laser flux density on 
the missile surface should be F =10 kW/cm* over a radius spot w=5cm, with an 
engagement dwell time of 3 seconds. Therefore, the average power delivered on the 
missile’s surface is 


P, =F mw” =(10kW/cm’) 25cm? ~800kW . (3.1) 


The major limitations of transfer of laser energy from the source to the target are 
the absorption and the scattering of the light traveling through the atmosphere. The 
following exponential law describes the power removed from the beam 


Phen, (3.2) 


where ©, is the extinction coefficient due to aerosols at sea level, P is the average power 
of the light beam that leaves the ship, and Pr is the average power of the light beam after 
traveling distance R through the atmosphere. By supposing a typical value for 
a, = 0.05km' and the missile is located at distance R = 6 km, then the average laser 


power that must leave the ship, according to the Equation (3.2), is 


0.05km7!6km 
( ) 


P =P,e“ =(800kW)e =IMW. (3.3) 


Equation (3.3) indicates that the required power at the ship is well within the 


capability of an FEL with an average power of 1 MW. 
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B. POWER SPECTRAL DENSITY (PSD) OF A SHORT PULSE FEL 

Free electron lasers, as seen earlier, produce their energy by passing a relativistic 
electron beam through a structured magnetic field. Electrons accelerated during passage 
through a magnetic field create the optical radiation. The proposed FEL weapon design 
uses micropulses of relativistic electrons. The electron micropulses have been accelerated 
by an accelerator with RF frequency ©. Thus, the FEL operates in a pulsed mode, 


producing an output light beam with pulse repetition frequency ©. 

The FEL used for this study has an average power P =1 MW, optical wavelength 
X =1.0um, pulse repetition frequency =750MHz, and pulse duration of 1 ps. A 
representation of the light pulses with a rectangular pulse shape in the time domain is 


shown in Figure 10, where 27; is the pulse duration, T =Q"' is the pulse period, and @ is 


the peak amplitude of the light electric field in a micropulse. 





fy =P Time 


Figure 10. Short Pulse Free Electron Laser. 


As seen in the previous section, the dwell time of a single engagement is 3 sec. 
Thus, for a single engagement, the beam consists approximately of 2 billion pulses. This 
is a very large number, and in order to simplify our analysis, a single engagement is 
considered to contain an infinite number of pulses. Using this approximation, the 


following equation describes the light signal in the time domain, 


n=co 


a(t)= ¥° dcos(2n f,t+0,)[ [(¢-27). (3.4) 


n=—0co 
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where f, is the light frequency that satisfies c=Af,, c is the speed of light, @, is the 


random phase of each optical micropulse, and II(r) is the function describing a 
rectangular pulse shape, 


0, |\>7, 


M=\) Wer, (3.5) 


In reality, each pulse would have a smooth bell shape and not the idealized 


rectangular shape described in Equation (3.4). In order to be more accurate in our 


analysis, the following cosine function c(t) is used to approximate the shape of each 


pulse, 
: H<7,(1-r) 
t}<i,U-r 
rail 
c(t)=41/2] 1+cos = ,  T(l-r)s||<7,(1+r) (3.6) 
r 
_ I'|>7, (1+r), 


where r is a roll-off factor which ranges between 0 and 1. When r equals zero, the 


function c(t) represents a rectangular pulse shape. Normalized pulse shapes for r=0, 


r=0.5, and r=1 are shown in Figure 11. 
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FEL Pulse Shapes 
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Figure 11. Plot of Equation (3.6) for r=0, r=0.5, r=1, and T;=0.5 ps. 


The Fourier transform of Equation (3.6) is given by [Ref. 9] 


c(f)= sin (2n fT, ) cost avai) = oT sine(2n fT) cos (2m ff) BD 
HS = A(Ayt) 1—(47fT, ) 
Replacing the function II(t) with c(t), Equation (3.6) becomes, 
a(t)= ¥° Gcos(2n ft+0,)c(t—nT), (3.8) 


which is truncated in order to be able to evaluate the power spectral density (PSD) for 
a(t) [Ref. 10] 


_ ¥. acos(2n ft+0,)C(t-n7), <5" 


ante: = (3.9) 


0, otherwise, 


where the subscript 7T' denotes the truncated version of a(t), equals to 


T'=2( N+1/2)T, and N is an integer. 
a2 


The average power P of the FEL pulse train is given by [Ref. 10] 


7 *(t)d ji =a rH =1MW (3.10) 
“Ty t rk Ae I) a ; 


where the bar represents the frequency average, a(t) is given by Equation (3.8), T has 
the value defined earlier, and A,.(f) is the Fourier transform of a,.(t). By combining 


Equations (3.8), (3.10), and solving the latter, the amplitude @ of each pulse is obtained, 


ee. 1 
T1-% 


In order to evaluate the PSD of the FEL short pulse train, the Fourier transform of 





(3.11) 


a 


Equation (3.9), A,.(f ), is first obtained 


Ar (f) =F| a, (N]=5 y [e*s (f-f,)te Ons (f+ Ff.) |*C(S) seen, 
or, 
An(f)=5 y je Pr rele( f— feeb OLE Ga) (3.12) 


where the notation x *y represents the convolution between x and y, C ( 3 ) is the Fourier 


transform of Equation (3.6) and is given by Equation (3.7). Then the PSD of a(t) 


becomes [Ref. 10] 


P,(f) = lim | 


T'— 0 


Aor 
ae (3.13) 


The substitution of Equation (3.12) into Equation (3.13) yields the final form of 


the single-sided PSD of a(t) 


B(f)=4 CBT: 


and by using Equation (3.7), 
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cos| 2nr( f —f,)T, | 
1-[4r(f-4,) TJ 


Am 2 
seen B14) 


P(f) sinc 2n ( f - f,)T; | 





where the amplitude @ is given by Equation (3.11). Now it is possible to calculate the 
average power of the FEL short pulse train a(t) in a finite frequency band [ f, fy| of the 


spectrum as, 


h 
P. 4, =| P(f)4- (3.15) 
fi 


Figure 12 represents the distribution of power in the pulse train as a function of 
frequency according to Equation (3.14), for a frequency range from f = 2.97x10'* Hz to 
f = 3.03x10'* Hz. The spectrum is plotted for three values of the roll-off factor: r=0 


(rectangular pulse shape), r =0.5 , and r =1 (cosine bell shaped pulse). 
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Figure 12. Power Spectral Density for an FEL Operating in Pulse Mode. 
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From Figure 12, it is obvious that the PSD for both the rectangular and the bell- 
shape pulses is centered at the fundamental frequency f,=3x10"Hz or A=lum. 


However, it can be seen that the rectangular pulse shape (r =0) spreads the power in the 
spectrum much more than the latter ones. The calculated average power that is contained 
in the fundamental lobe, which corresponds to the spectrum [2.99x10'4 Hz, 3.01x10" 
Hz], is 90% of the total average power for the rectangular case, while for the bell shaped 
cases, it is more than 99.9%. Thus, the cosine bell-shaped pulse decreases the sidelobes 
of the power distribution in the frequency domain and concentrates most of the power in 
the main lobe close to the fundamental frequency, which is desirable for laser beam 


propagation through narrow transmission windows in the atmosphere. 


In conclusion, the FEL operation in pulse mode leads to power spread in the 
frequency domain, depending on the FEL pulse shape in the time domain. The 
propagation through the atmosphere of a broad-spectrum laser pulse may have several 
consequences: 


e Presence of power in frequency components where the absorption 
coefficient is high. This will increase the total absorbed power and also 
may cause nonlinear effects such as thermal blooming. 


° Presence of a sufficient amount of power in the visible region, where the 
scattering process may make the laser beam visible. In most cases, in a 
real operational environment, this will be undesirable for a weapon 
system. 


e Presence of a sufficient amount of power in the region of the spectrum that 
is not considered eye safe. This will compromise the personnel safety on 
the ship, especially those who will have eye contact with the beam director 
during the engagement. 


e Finally, the presence of a sufficient amount of power in the microwave 
regime that may lead to a soft kill mechanism of an incoming missile 
threat by damaging the electronics the missile. 


In the present study, the possible problem of thermal blooming for the proposed FEL 


weapon design will be examined. 
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Cc. THERMAL BLOOMING 
The issue of whether the propagation of a MW- level FEL laser beam through the 


atmosphere causes thermal blooming will be studied next. An FEL operating in pulse 
mode, as previously seen, leads to power spread in the frequency domain. The worst-case 
scenario is that from a square optical pulse shape. Figure 13 shows the normalized power 


spectral density of such a pulse train for the spectrum [10'* Hz, 3x10'° Hz]. 


Normalized Power Spectral Density 
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Figure 13. Normalized power spectral density of a square optical pulse shape. 


The absorption coefficient profile versus frequency for a moderately clear 


maritime atmosphere with visibility of 23 km is shown in Figure 14. 
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Absorption Coefficient Profile 
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Figure 14. Atmospheric Absorption Coefficient Profile Versus Frequency for 


Moderately Clear Atmosphere (23 km Visibility).(After [Ref. 7]). 


The power that is contained in the spectrum of Figure 13 equals to 999.85 kW. 
The propagation of that amount of power through the atmosphere will lead to a specific 
amount of absorbed power. Consider one liter of air that has about the same cross section 
as the laser beam. Then, based on Equation (3.2) for R=10 cm and the absorption 
coefficient profile of Figure 13, the absorbed power by | liter of air is calculated to be 


0.255 W, or an absorbed energy rate AE/At =0.255 Joule/sec. For a CW FEL (no pulse) 


operating at the same frequency and at the same atmospheric conditions, the absorbed 
power would be 0.25 W. From these results, it can be seen that the absorbed power 
increases only by a small amount (0.005 W) due to spectral spread of the short optical 


pulse. 


The specific heat of air C, is given by 


AE J 
C,=—— = 1 ; 
VAT IC 





, (3.16) 


a7 


where V is the air volume and AT’ is the temperature rise of that volume. By considering 
a dwell time of 1 sec and solving the last equation for AT the rate of temperature rise of 


the air volume equals to is 


°C 
— ———~ = 0.25—. (3.17) 
At VC, S 
LJ} 1 





AT AE/Ar 0.25 Y 
J 


0 


LC 


Taking into account the ship’s movement, the incoming missile’s velocity, and 
the local wind direction, a stagnate air column (stagnation zone) can exist somewhere 
along the laser beam path. Statistically, for a dwell time of 1 sec, the size of a stagnation 


zone may vary from Az =1m to 100 m [Ref. 7]. 


Any small change in the temperature of the stagnate air column causes a small 


change in the index of refraction of the air. The index of refraction n is given by 


n=—, (3.18) 
v 


Cc 


where c is the speed of light in vacuum and v, is the speed of light in air. Then, the small 


change in the index of refraction will be 





An =-——*+ = -—+, (3.19) 


since v. =c. By relating An to the temperature change AT , [Ref. 7] 


An=-10°AT(°C). (3.20) 


From the last two equations, the change in the speed of light is 


Av, =-c An =-10°cAT (°C). (3.21) 


This leads to a lens effect. The inner core of the stagnation zone will be hot and 
the light travels at speed v, +Av,, while the outer part of the zone is less warm and the 
light travels at speed v. That difference in light speed leads to a difference in path length, 
where the light at the inner part of the beam gets ahead of the light at the edge of the 


beam by 
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As = Av, 92 =—Andz=10°AT (°C) Az, (3.22) 
c 


after traveling through a stagnation zone of length Az. The racing of the light traveling 


through a stagnation zone is shown in Figure 15. 
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Figure 15. Light Traveling through a Stagnation Zone. 


The effect of the stagnation zone is equivalent to a negative lens. Thus, thermal 
blooming is a result of the atmospheric heating due to absorption, which creates a 


negative lens in the atmosphere and spreads the beam, as shown in Figure 16. 
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Figure 16. Laser Beam Passing through a Negative Lens.(After [Ref. 4]). 





The lens effect is measured by comparing the path difference As to the optical 


wavelengtha . Then, from Equation (3.22) 
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As Az 
—=10°AT (°C). 3.23 
5 ( ls (3.23) 


Considering a dwell time of lsec, A =lum, Az=10m, then AT =0.25°C and Equation 


(3.23) becomes 


AS 10 (0.25) 
n 


a =2.5. (3.24) 

If a stagnation zone of length Az=50m is considered, then As/A =12.5 
wavelengths. Thus, it can be concluded that the blooming effect of a MW level FEL at 
A =1um is not strong for a stagnation zone size less than 10 m, but it gets quite strong 
for greater values of Az. Of course, at this point, it is necessary to mention also that our 
approach is only based on the atmospheric absorption profile for a moderately clear 
atmosphere. The results may vary for different atmospheric conditions. 

Another significant parameter is the position of the stagnation zone along the 
beam path. Since the stagnation zone tends to spread the hser beam core, even a strong 
blooming effect close to the target may be preferable to a small stagnation zone that is 
positioned close to the beam director. 

Also, there are more detailed transmission windows in the atmosphere that are not 
included in Figure 14. The result for a window with higher transmission would decrease 
absorption and blooming. 

D. THE EFFECT OF THERMAL BLOOMING ON THE LASER BEAM 

Beginning from Maxwell’s wave equation, the dimensionless parabolic wave 
equation in free space (Equation 2.50) can be obtained which describes the diffraction of 
the laser beam from the beam director to the target. One solution of Equation (2.51) is a 


Gaussian laser beam described by Equations (2.53) and (2.54). 

As seen previously, the atmospheric absorption in a stagnation zone causes 
thermal blooming. Considering that the beam during the propagation to the target hits a 
stagnation zone at a position defined by the dimensionless timet,, then the phase 


distortion caused at the beam by the thermal blooming is described by [Ref. 7] 


Aout (x, y,t,) = nes Pg (x, yt. ) 2 (3.25) 
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where a,,, is the field entering the zone, a,,, is the field after the zone, and d9, is the 
strength of the zone lens. The phase change factor of the last equation can be related to 


the lens effect measured by the path difference As over the optical wavelength 1, 


2 As 
lacy.7 | dp, =. (3.26) 
Using a simulation code, it was possible to solve Equation (2.50) numerically for 
the Gaussian laser beam described by Equations (2.53) and (2.54). The program follows 
the evolution of the laser field over many small time steps dt ina transverse window 
consisting of nx bin elements in both the x and y directions. The simulation shows the 


effect of the diffraction on a laser beam that is focused on the target with or without 
thermal blooming. 

The results of the simulation are plotted in Figures 17 and 18. The color intensity 
of each plot represents the field intensity strength at each bin as shown at the bottom of 
each plot. Each figure contains two plots. The plot on the left is a side view of a slice 


through the laser beam as the feld a(x,nx/2,t) evolves in dimensionless time from 
t =0 to 1. The plot on the right side represents the spot of the laser beam at the target, 


Ja(x.y.t 7 1). 


Both plots in Figure 17 indicate a focusing beam due to the diffraction in the 
absence of thermal blooming, while Figure 18 shows a distortion produced at the laser 
beam by the stagnation zone of strength dd, =10 at the position Tt, =0.7. Since the 
dimensionless optical field amplitude was chosen to be unity, the stagnation zone has a 
length of about 50 m. The stagnation zone creates a “hole” in the laser beam, which is 


obvious in both the side and the front view plots. 
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Figure 17. Laser Beam Focused on the Target Due to Diffraction. 
Side View for y = nx/2 Front View at r=1 
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Figure 18. Laser Beam Focused on the Target Due to Diffraction While Hitting a 


Stagnation Zone at tT. =0.7. 
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E. CONCLUSIONS 

In this Chapter, the spectrum of the laser beam output from the proposed FEL 
weapon design is obtained. The proposed FEL operates in a pulse mode, producing 
optical micropulses 1 ps long at a pulse repetition frequency of 02=750 MHz. The 
spectrum may cover many absorption lines. The absorbed laser energy heats the air and 
may lead to nomlinear effects such as thermal blooming. It is found that thermal 
blooming is not an issue for a moderate clear atmosphere when the stagnation zone size 


remains less than 10 m. 


The next Chapter presents the simulation results of the proposed FEL weapon 
system. The effects of the electron beam focusing and electron beam off-axis shifting are 


studied. 
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IV. SIMULATIONS OF THE 1 MW WEAPON DESIGN FEL WITH 
STRONGLY FOCUSED ELECTRON AND OPTICAL BEAMS 


A. SHIPBOARD FEL COMBAT SYSTEM 
A Naval FEL application is defense from anttship missiles at ranges around 
5 km. As seen previously, the destruction of a missile at that range requires a laser beam 


of about 1 MW power at around 1 um wavelength for an engagement that lasts a few 


seconds. 


The naval application also requires a compact configuration with the laser 
components packaged into a rectangular module placed onboard a ship. The proposed 
FEL MW design has an estimated size of 12 m x 4m x 2 m and is shown in Figure 19. 
The figure represents a top and a side view of the module. The accelerator runs along one 
side of the box while the optical resonator lies on the opposite side, leaving the support 
components in the area between these two. The whole system is estimated to weigh about 
40 tons with an estimated cost of about $60 M [Ref. 11]. A possible arrangement of the 
shipboard FEL system, as presented during the 2001 FEL workshop in Newport News 
VA, is shown in Figure 20 [Ref. 1]. The modular support systems, such as the RF 
klystrons (red boxes), are located separately but near the other components. This layout 
shows the system in a horizontal orientation near the keel penetrating a single water-tight 


bulkhead, which is an undesirable but tolerable feature. 
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FREE ELECTRON LASER 





<< | 12 
Figure 19. Proposed FEL MW Design. (From [Ref. 11]). 








Figure 20. Conceptual FEL Weapon System Ship Arrangement (red boxes: RF 
klystrons, RF power components - blue boxes: refrigerator, accelerator, wiggler — 
purple boxes: bending magnets). (From [Ref. 1]). 


The requirement of high output power from a small compact FEL onboard a ship 
raises some challenging issues, such as the optical intensity on the mirrors and 


vibrationally induced misalignments of the electron beam and the cavity mirrors [Ref. 1]. 
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1. Short Rayleigh Length 

The large intracavity power and a small optical beam size lead to high intensity on 
the cavity mirrors, making them vulnerable to damage. This problem could be solved by 
inventing more robust mirrors, by lengthening the cavity n order to increase the mode 
area at the mirrors by diffraction, or by using a shorter Rayleigh length [Ref. 12]. In the 
current design, however, the resonator length is limited to S = 12 m. Additionally, it is 
estimated that cooled sapphire mirrors with a transmissive coating can handle an intensity 
of ~200-300 kW/cnY [Ref. 11]. Thus, the use of a short Rayleigh length resonator is 
required to increase the spot size on the mirrors, leaving a very small optical mode waist 


for the intense FEL interaction. 

The requirement of a short Rayleigh length resonator for a high power FEL can be 
derived as follows. A typical FEL has an undulator length of L=3m and a Rayleigh 
length of Z, =~ L/3=1m. For 4 =1 um, the optical mode waist according to Equation 


(2.51), Z,A=nW? 


o°? 


would be about 0.6 mm. Then, from Equation (2.52) for 
z=S/2=6m, the laser spot size at the mirrors would still be a few mm (W ~3.6 mm). 
Assuming 25% output coupling, 1 MW output power requires a power stored inside the 
cavity of 4 MW. The intensity on the mirrors would then be about 10 MW/cnr, which, of 
course, will lead to mirror damage. By reducing the Rayleigh length to Z, =1.8 cm, the 
spot size of the beam at the mirrors is 2.5 cm, and the intensity is reduced to 210 kW/cnr. 
This is within the tolerance of the mirrors. However, for Z,=1.8 cm, the optical mode 


waist is very narrow in the center of the cavity, W, =0.08 mm. 


oO 


A short Rayleigh length implies a very small optical mode waist and a beam 
which expands rapidly towards the end of the undulator. In the case where the electron 
beam radius is greater than the mode waist radius, r, >W,, many electrons will not 
overlap with the intense optical field in the interaction region near the cavity center. 
Thus, a narrow electron beam is desirable. Additionally, the undulator length must be 
short in order to avoid damaging the undulator magnets by the intense rapidly-expanding 


laser beam. 
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2. High Power FEL Parameters 
The goal is to design an FEL with 1 MW output power, using a short Rayleigh 
length optical cavity. Therefore, the optical resonator could have the characteristics 


shown in the table below. 


OPTICAL RESONATOR PARAMETERS 


Resonator Length S=12m 
Optical Wavelength A=lyum 


Quality Factor Q,=3 (28% power transmission per 


pass) 


Rayleigh Length Z, =1.8 cm 


Optical Resonator Parameters fora MW Level FEL.(From [Ref. 13]). 





As previously seen, the short Rayleigh length limits the undulator length. In this 
case, the linearly polarized undulator has a length of only 1/20" of the resonator length. 
In this manner, fewer undulator periods are used, which enhances the extraction 
efficiency according to Equation (2.33). The magnetic undulator has the characteristics 


shown below. 


MAGNETIC UNDULATOR PARAMETERS 


Undulator rms Parameter 


Magnetic Undulator Parameters for a MW Level FEL. (From [Ref. 13]). 
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Finally, the MW FEL design calls for an electron beam with the characteristics 


shown below. 


ELECTRON BEAM PARAMETERS 


Pulse Length 


Extraction Efficiency 


Electron Beam Parameters for a MW level FEL. (From [Ref. 13]). 





B. BETATRON MOTION 
The FEL weapon design uses a linearly polarized undulator. Thus, it is necessary 
to calculate the electron motion in this kind of undulator. The undulator magnetic field 


can be written as 


B, = B[0,sin (k,z)cosh(k, y),cos(k,z) sinh (k, y) ], (4.1) 


where k, =2n/X, is once more the undulator wavenumber. The magnetic field near the 
undulator z-axis (y ~0) is 


B, = Bsin(k, 2) 3. (4.2) 


By squaring and averaging over an undulator period to remove the fast electron 


wiggling motion Equation (4.1) we have 
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[ cosh? (k,y)+sinh? (k,y) | 
a5 s+ sinh? =) (4.3) 
~ B? E +(k,y)" tof 


Thus, the rms value of the magnetic field is given by 


: ~ A (1+ (k,y) +). (4.4) 


V2 


From the last equation, it is obvious that the field increases off-axis with y, 


Bins = B. 








focusing the electrons back towards the undulator axis. Near the undulator axis, 


k,y = 0, the undulator field is approximately B,,. = B/ oe 


The total force acting on the electron comes from the Lorentz force Equation 


(2.13), with the magnetic field given by Equation (4.1). Assuming no light is present 











(y =0,£, =0), 
— ee 5 [> cos(k,z)sinh(k,y)—z sin(k,z)cosh(k,y) |, 
ym,c 
B, = a icos(k,z)sinh(k,y), (4.5) 
; CB pcs 
6. =- > xsin(k,z)cosh(k,y). 
ym,.c 


The first equation can be solved exactly by integration. Assuming perfect electron 
injection near the undulator axis in order for the integration constant to be zero, it 
becomes 


B 
B, =- Bae cos(k,z)cosh(k,y) 


eae V2 eB AX, 


=e cos(k,z) cosh(k,y) (4.6) 


2K 


=— Cos (k,z)cosh(k,y), 
Y 








50 


where K is the undulator parameter, given by Equation (2.3). Combining the last two 


equations we have 
. cK°k : 
Biro eam eEY), 


2 
6. = ae sin (2k,z)cosh’ (k,y). 
Y 


(4.7) 


For an electron near the undulator axis, (k, y <<1) , the electron motion over 
many undulator periods can be determined. Averaging Equation (4.7) over many 1,, 
gives 6. =0 or B. = B,, and 6, = -(cK “k /y’) y. Therefore, the averaged equation of 


motion in the y direction is 


2721712 2 
jones ki sinh(2k,y) =- cKk, y. (4.8) 
2y Y 








Again introducing the dimensionless time t=ct/L and the dimensionless 


derivatives (...) = d(...)/dt , the last equation becomes 


2 
_ Kk,L 
p=-{ Be) y. (4.9) 


From Equation (4.9), it is obvious that the electron motion in the y direction is that of a 


simple harmonic oscillator, with a dimensionless betatron frequency 


_ Kk,L _ 2nNK 
oy Y 








(4.10) 


For the parameters of the weapon design @, = 2m /8, the electron executes only 
1/8 of a betatron oscillation along the undulator axis. Thus, the electron undergoes a slow 


oscillation with frequency @, in the y direction along the undulator axis. 


A general solution of Equation (4.9) has the form 





v(t) =), cos(oyt) + sin(@,t ), (4.11) 
B 
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where y, is the initial (t =O) offset y-coordinate of the electron entering the undulator, 
8, =y(0) / L is the initial angle of the electron velocity with respect the undulator axis on 


the yz plane, and y(0) is the initial electron y-component velocity. 


The betatron motion of an electron with initial coordinates y, and 8, in the y 


direction was studied. Since the field is independent of the x coordinate, the trajectory can 
be displaced in the x direction without any change in B. or the resonance condition. In 


an experiment, any displacement of the electron trajectory in the x direction is confined 
by the extent of the optical mode, or by using additional focusing magnetic fields. 


However, an angle 8, does affect the resonance condition according to Equation (2.31) 
by changing B. for a given y. Since an electron’s phase velocity is given by 


veL [ (« +k,)B. -k] , it can be shown that an electron injected with a displacement y, 


and an angle 8, =,/0; +®; reduces the phase velocity by [Ref. 5] 





ment 


Av= 5 
1+K 


K*k?y +702). (4.12) 


ovo 


In reality, the beam consists of a very large number of electrons that enter the 


undulator with a random distribution of displacements and angles with mean values ¥: 


and 6, . The product of these two mean values defines the beam emittance. In most cases, 
the contributions of these two values are matched using focusing magnetic fields, so that 
the beam does not expand too much as it moves along the undulator and the phase 
velocity spread is minimized [Ref. 14]. 
C. THREE - DIMENSIONAL SIMULATIONS 

It is important that the shipboard FEL combat system be isolated from shock or 
vibrations from the ship’s deck and hull. The vibrations may lead to an electron beam 
misalignment, either an offset or a tilt, or to a cavity mirror misalignment. Given a small 


optical mode waist radius (W, =0.08 mm), a slight misalignment of the electron beam 


7 


could conceivably reduce the overlap between the electrons and the intense optical field 


in the interaction region. 


A given accelerator has a fixed normalized beam emittance, 


€ =yr8o, (4.13) 


where 0, is the electron beam rms angular spread, and 7, is the rms electron beam radius. 


In most cases, re and ®, are matched so that the beam does not expand too much as it 
moves along the undulator [Ref. 13]. For a typical Rayleigh length cavity, a matched 
electron beam optimizes overlap with the optical mode. In a short Rayleigh length cavity, 


the electron beam could be focused with additional external magnets along the undulator, 


decreasing 7. and increasing 9.. The electron beam focusing optimizes the overlap with 
the optical mode, which has a very small waist in the interaction region, and thus 
enhances the gain and efficiency. However, the induced energy spread in the electron 
beam increases and is given [Ref. 5] 
AE _Ay_ Av 
E y 4iN’ 





(4.14) 


where Av is the induced electron phase velocity spread due to the FEL interaction. 


In the simulations presented here, the effects of the electron beam off-axis shift in 
the y direction for a high power FEL with the parameters mentioned in the Tables 1 
through 3 are studied. Figure 21 indicates the geometry of the problem. The figure shows 
the optical mode (blue lines) of a short Rayleigh length undulator inside the optical 
cavity. The electron beam (red line) is shifted from the undulator axis (green line) by a 
distance y,. The experimental design tolerances assume the beam will be aligned to 


within 0.01 mm of the undulator axis [Ref. 13]. 


Additionally, the electron beam focusing is studied. Figure 22 shows that the 
electron beam (red line) is focused to overlap the narrow and intense optical mode (blue 
line) in the center of the undulator. During the focusing process, the electron beam waist 
radius r. is reduced, while the angular spread of the beam is increased in order to keep the 


emittance constant. 


a1) 


E-beam 


Optical Mode 


Undulator 





Figure 21. Electron Beam Shift. 
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Figure 22. Electron Beam Focusing. 


1, Simulation Methods 

In a real FEL experiment, an optical wavefront with finite transverse dimensions 
will have multiple transverse modes determined by the interaction with the electron beam 
and the mirrors of the optical cavity. The three-dimensional simulations in x, y, and ¢ 
include the effects of diffraction, optical mode distortion, mirror transmission and edge 
losses. The undulator is oriented along the z-axis, with the magnets normal to the y-axis. 
The simulations use dimensionless coordinates. The longitudinal lengths in the z direction 


are normalized to the undulator length L, the transverse lengths in x and y directions are 


divided by the characteristic mode radius (Li/n)”, and the angles are divided by 


(A/m aa . Thus, the dimensionless parameters for the high power FEL presented here are 
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DIMENSIONLESS PARAMETERS 


Electron Beam Radius ; ; a i LA/n i 03 


Rayleigh Length z, =Z,/L= 90.03 


Electron Beam Parameters for a MW level FEL. (From [Ref. 13]) 





The electrons are given an initial spread in positions (@,y) and angles (0,8, ) 


determined by the beam radius, beam emittance and focusing. They may also have an 
offset in position or angle, to study the effects of beam misalignment. They evolve in 
longitudinal phase space according to the FEL pendulum equation (2.28). In addition, 
they undergo betatron oscillations in the yz-plane. From Equation (4.11), which describes 
the betatron oscillation for a single electron 

0, +A0, 


y(t) =(y, +Ay)cos| @, (« -t, 5 (t-t, )]; (4.15) 


where y, is the initial electron beam offset, 8, is the initial electron beam tilt, Ay is a 


random spread in electron position due to emittance and the beam radius r., and A@, is 


an electron random angular spread due to emittance also. The dimensionless time 


t =ct/L corresponds to the electron’s position along the undulator axis, and T, is the 


position where the electron beam is focused to its minimum size. For the present study, 


the electron beam is shifted and focused at the undulator center with no tilt, y, #0, 
8, =0, and t, =0.5. The electron phase velocity will be reduced due to the beam 
misalignment according to Equation (4.12). 


The optical wavefront evolution is described by the parabolic wave equation, 


starting with an initial Gaussian profile and evolving over many passes until a steady- 


a2 


state mode is obtained. Then the single-pass extraction efficiency is calculated using 


Equation (2.33), 


=-er——— 4.16 


where (Av) is the shift in average electron phase velocity due to the FEL interaction. 

2 Three-Dimensional Simulation Output Format 

The output from a three-dimensional simulation of the proposed MW level FEL, 
using a short Rayleigh length of z,=0.03, is shown in Figure 23. The dimensionless 
parameters for this experiment appear in the upper-right window. The electron beam has 
a current density j = 210, a radius in both x and y directions of 0, , =0.3, and for this 
experiment, there is no electron beam shift or tilt (so y, =0, 0 —- 0). At each pass, the 
electrons are injected with an initial phase velocity v, =10.66 and zero initial Gaussian 


spread in phase velocity (6, =0), which optimizes the extraction efficiency (N =2.5% ). 


FEL WAVEFRONTS 


j=210 o =0.3 


Vv =10.660.=0.3 
@,=0.7 oc, =0.04 
By 
T,=0.5 Go, =0.04 
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° vy 
86 =0 o_=0 
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z =0.03 r,=0.0 
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m c 
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Figure 23. Three-Dimensional Simulation Output for the Proposed MW Level FEL 
Using a Short Rayleigh Length with a perfectly-aligned Electron Beam. 
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The electron beam is focused at the undulator center (t, =0.5) and has an 


angular spread 6, =0.04, which is determined by the matching requirement 


Oo =0;0;, [Ref. 14]. The dimensionless mirror radius is r, =13.2, which allows 1% 


(e=0.01) edge losses per pass, and the dimensionless radius of curvature is r, =1.5, 


determining the Rayleigh length [Ref. 14]. The color intensity scale, shown on the right, 
is used for the optical field amplitude intensity plots. The scale runs from dark blue, 
corresponding to zero amplitude, to light blue, corresponding to maximum amplitude. 


The white contour line corresponds to a field amplitude at 5% of the maximum value. 


The upper portion of the figure shows two intensity plots of the optical field. At 


the upper left is the side view of a slice through the optical wavefront as the field 


la (y .n) evolves over n=140 passes. The one in the upper center represents the spot of 


the final wavefront at the undulator exit le(x.¥) after n=140 passes. In the same 


wavefront plot, the red dots show the position of sample electrons, indicating any offset 


or misalignment of the beam. The center plot presents a side view of a slice through the 


optical wavefront a(x.) during the final pass. The red dots correspond to sample 


electrons as they move from the beginning of the undulator (t =0) to the end (t =1). 
The optical wave amplitude profile, corresponding to the power distribution at the 
mitrors, is shown with the yellow contour lines. In this simulation, the mirror separation 
is only three times the undulator length, for numerical convenience. 

The simulation starts with a weak optical field amplitude of a, =5 and after 
about 100 passes saturates to its final value of la| = 69. This is also shown in the lower 
right plot of the optical power P(n) and gain G(n) evolution, where the oscillator 
saturates after about 100 passes. The lower left plot is the evolution of the electron 
spectrum f (v,n). Here, it is shown that the initial electron phase velocity is v, =10.66. 


It is also shown that the energy spread of the electron beam increase due to the FEL 
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interaction until the laser saturates. The last plot in the lower center shows the final 
electron phase-space evolution where evidence of bunching is shown. 
3 Simulation Results 
a. Electron Beam Shift 


The simulations for the MW level FEL design are run using a short 
Rayleigh length for various electron beam shifts y,. The normalized electron beam shift 
was varied from y,=0 to y, =1 in steps of 0.1. As mentioned previously in this chapter, 
the experimental design tolerance assumes that the beam will be aligned to within 0.01 
mm of the undulator axis. The normalized beam shift that corresponds to 0.01 mm is 
y, =0.023 . At each value of y,, the initial phase velocity v, was also varied, in order to 
determine the peak efficiency 1. The efficiency is calculated according to Equation 


(4.16) after steady-state is reached. 


The steady-state FEL efficiency 1 versus the initial phase velocity for 


four values of the normalized electron beam shift ( y, =0,0.2,0.4,0.6) is plotted in Figure 
24. 


58 


Efficiency, sz (%) 





-10 -5 0 5 10 15 20 
Initial Phase Velocity, v. 


Figure 24. Single-Pass Extraction Efficiency n Versus the Initial Phase Velocity v, 
for Four Values of the Normalized Electron Beam Shift y,. 


For each value of the electron beam shift, the efficiency steadily increases 
as the phase velocity increases, with a sharp drop-off just after the peak. The peak 
corresponds to the value of v, where the FEL gain drops below threshold. It can be seen 
that with no ekctron beam shift (y, =0), the FEL has the highest peak efficiency of 
1 =2.5%, which corresponds to initial phase velocity v, =10.66. This was expected to 


be the peak, since for zero shift, the electron beam is almost perfectly aligned with the 
optical mode. The peak efficiency decreases as the electron beam is shifted from the 
undulator axis. Additionally, the electron beam shift changes the initial phase velocity for 


which the laser reaches the peak efficiency. However, the phase velocity varies with the 
optical wavelength as Av=2nN(AA/A). Thus, an electron beam shift can change the 


optical wavelength. 


ao 


Figure 25 shows the output of a single simulation with the electron beam 
shifted off-axis by y, =0.6. The FEL interaction takes place close to the undulator center 


as the electron moves toward the end of the undulator. The electron beam is not aligned 


with the intense optical mode in the interaction region due to the beam shift. This causes 


the optical mode to tilt, which is obvious in the la(y.t ) plot. Nevertheless, the optical 


power P(n) reaches steady-state, but the extraction efficiency is reduced from 2.5% to 


1.8%. The phase velocity is reduced according to Equation (4.12) by Av= 0, y, =(0.2. 
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Figure 25. Three-Dimensional Simulation Output for the Proposed MW Level FEL 
Using a Short Rayleigh Length with an Electron Beam Shift y, =0.6. 
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Figure 26 summarizes the results of the simulations. It shows the peak 
single-pass extraction efficiency versus the normalized electron beam shift varied from 


y, =0 to y, =1 insteps of 0.1. 
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Figure 26. Peak Single-Pass Extraction Efficiency 1] Versus the Normalized Electron 


Beam Shifts y,. 


In this figure, the normalized shift that corresponds to 0.1 mm is indicated 
for reference by an arrow on the horizontal axis at y, =0.23. In the same way, the 
normalized experimental design tolerance that corresponds to 0.01 mm beam shift is 
indicated at y, =0.023 . The required efficiency to achieve the 1 MW goal is indicated by 
a dashed line at rj =0.7 %. Here, it is once more clear that the peak efficiency steadily 
decreases, as the beam is further offset from the undulator axis, but remains above the 


goal for y, <0.75(=0.3mm). This value is well beyond the design tolerance of 
y, = 0.023(0.01mm). 


Additionally, as was presented by the NPS FEL group at the 25" 


International FEL Conference of 2002 in Argonne, IL, the electron beam tilt (6, #0) 
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with no offset shift (y, =0), gives similar results [Ref. 13]. The efficiency remains above 
the 0.7% extraction efficiency goal for 0, < 8(6mrad) when the beam is tilted about the 
undulator center (t, = 0.5), and for @, < 1.5(1mrad) when the beam is tilted about the 
undulator entrance (t, = 0). Both values of the electron beam tilt are well beyond the 
design tolerance of 20 prad , which corresponds to normalized 8, =0.03. Thus, it can be 


concluded that the shipboard MW level FEL design is unaffected by any electron beam 
misalignment (an offset or a tilt) due to vibrations, if the system is isolated within the 
previously mentioned design tolerances. 

b. Electron Beam Focusing 

All the results presented so far have used a matched electron beam, 
without using any focusing to enhance the overlap between the electron beam and the 


optical mode. The beam can be focused, by varying the electron beam waist radius 7,, 
while keeping the total current and emittance fixed. In the simulations, the normalized 


electron beam radius 0, = 1, /(LA/x ” is used. 


Figure 27 shows the simulation results for the peak efficiency 1] versus 
the normalized electron beam radius, which is varied from o, =0.3 to 0.15 in steps of 


0.05. The normalized beam radius that corresponds to 0.1 mm is indicated by an arrow on 


the horizontal axis at 0, =0.23. The largest value of 6 , =0.3 corresponds to a matched 


beam for the parameters of this FEL design. As 6, is reduced by focusing the beam, 


the peak efficiency increases from 1} =2.5% to 4 %, which is far beyond the goal value 
of 0.7 %. This indicates that it may be possible to reduce the average current and still 
obtain the goal of 1 MW output power. In the same figure, the induced energy spread is 
also shown, which is given by Equation (4.14). The energy spread needs to be kept below 
the maximum allowable limit for safe electron beam recirculation of about 15 % [Ref. 
15]. The results indicate that the induced energy spread increases from 11% to 14% as the 


beam is focused from 0, =0.3 to 0.15. 
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Figure 27. Peak Single-Pass Extraction Efficiency rj and Induced Energy Spread 


AE/E Versus the Normalized Electron Beam Shifts y,. 


Figure 28 shows the output of a single simulation where the electron beam 
is focused to 0, = 0.15. Since the electron beam radius has been reduced by a factor of 2, 


the electron density p is increased by a factor of 4. The dimensionless beam current 
density (jp) is likewise increased by a factor of 4 to j =840. The electron beam 
angular spread is also increased according to the requirement Oo = O59; =0.16 to 
keep the emittance constant. Beginning with an initial velocity v,=12.8 and initial 
optical field amplitude of a, =20, the optical power P(n) reaches a steady-state after 


= 50 passes. 
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Figure 28. Three-Dimensional Simulation Output for the Proposed MW level FEL 
Using a Short Rayleigh Length and Focused Electron Beam. 
The initial phase velocity increased from v,=12.8 to 13 when the 
simulation was run again as presented in Figure 29. This change represents a very small 


change (0.1%) in optical wavelength. This time the optical mode la(y.t) , the optical 





power P(n), and the extraction efficiency oscillate over many passes, without reaching a 


steady-state condition. The complicated pattern of the la( y,t ) plot shows the presence 


of multiple transverse optical modes in the cavity. The laser has not reached a steady- 


state after even 600 passes, but expresses an extended oscillatory behavior. 
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Figure 29. Three-Dimensional Simulation Output for the Proposed MW Level FEL 
Using a Short Rayleigh Length and Focused Electron with Extended Oscillations. 


D. CONCLUSIONS 

The proposed FEL weapon design uses a short Rayleigh length undulator to 
increase the optical beam spot size at the mirrors. But a short Rayleigh length implies a 
small optical mode waist in the center of the cavity. A slight misalignment of the electron 
beam reduces the overlap between the electrons and the optical mode in the interaction 
region, resulting in less gain and efficiency. Using three-dimensional simulations, it is 
found that the shipboard MW level FEL design is unaffected by any electron beam 
misalignment (an offset or a tilt) due to vibrations, if the system is isolated within the 
expected design tolerances. It is also found that that the electron beam focusing enhances 


the FEL extraction efficiency. 


Next chapter presents the simulation results of the upgraded to 100 kW Thomas 


Jefferson National Accelerator Facility (TJNAF) FEL using a step-tapered undulator. 
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V. SIMULATIONS OF THE 100 KW TJNAF FEL USING A STEP- 
TAPERED UNDULATOR 


A. INTRODUCTION 

The Thomas Jefferson National Accelerator Facility (TJNAF) has developed and 
demonstrated the most powerful FEL in the world, with an average output power of 2 kW 
that is tunable in the infrared. This output power is not considered enough for a military 
DEW weapon system. But a major upgrade of the TJNAF FEL system to 10 kW is 
underway, and further extension of the system to 100 kW has already been proposed. 


This begins to open a credible path to the desired MW level FEL [Ref.1]. 


A schematic representation of the TJNAF FEL was shown in Figure 5 of Chapter 
II. The system recirculates the electron beam to recover the excess energy of the electron 
beam. It extracts the maximum energy from the electron beam in a single pass through 
the undulator while inducing the minimum amount of energy spread for beam 


recirculation (energy spread Ay/y <15%). 


At the Naval Postgraduate School, in collaboration with TJNAF, we have 
explored the use of a step-tapered undulator. This undulator design alters undulator field 
in order to alter the resonance condition halfway through the undulator to improve the 
performance of the proposed 100 kW TJNAF FEL. The use of short electron pulses 
complicates the desired interaction. The results of the experiment were presented during 
the 23" International FEL Conference of August 2001 in Darmstadt, Germany [Ref. 16]. 
B. STEP-TAPERED UNDULATOR 

Unlike the conventional undulator, which has a periodic field and wavelength, the 
step-tapered undulator abruptly changes the value of the field and therefore the undulator 
parameter K halfway through the undulator. This could be done by having two separate 


undulator segments, each with different value of the magnetic field and K. 


A step-tapered undulator results in a modified pendulum equation [Ref. 5, 14-16] 


oe oft 5 ptacost +o), (5.1) 
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where C =(k+k,)z—@t is the electron phase, k is the optical wavenumber, k, = 21/2, 
is the undulator wavenumber, @ is the optical frequency, v is the electron phase 
velocity, @ is the optical phase, a is the dimensionless optical field amplitude, @(z) is 


the step function, 


0, z<0 


: 52 
1, z20 ef) 


(2)={ 


and A is the phase acceleration caused by the step taper AK/K at t =z/NA, =0.5 along 


4n NK* \( AK 
fee |e, (5.3) 
1+K K 
Thus the stepped field is decreased when A> 0, or increased when A<O. 


& 100 KW TJNAF FEL PARMETERS 
The TJNAF FEL can be upgraded to operate at increased output power of 100 kW 


the undulator and is equal to 


by increasing the kinetic energy of the electron beam from the accelerator to 


K, =210MeV . The electron beam consists of pulses with length 1, =0.lmm, pulse 
repetition frequency Q =750MHz, and peak current 1=270A. The average electron 
beam power would then be 


pete amiw. (5.4) 
Cc 





So an output optical power of 100 kW requires an extraction efficiency of n ~ 0.7%. The 
undulator period is 2,=8cm over N=36 periods with an rms undulator parameter 
K=1.7. 

The electron beam is described by dimensionless peak current j=5 and 
dimensionless electron pulse length o, =/,/NA =3. The resonator quality factor is 
Q, = 4.2, which corresponds to cavity losses per each pass of 1/0, = 24%. 

D. STEADY STATE POWER 


To study this FEL, we run the simulations starting with a weak optical field 


(a<m) for many passes until the power reaches a steady state condition. The step-taper 
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is varied from A=—2n to A=2n in steps of 1. At each value of A, the desynchronism 
d =-—2AS/NxX is also varied, in order to determine the peak efficiency (Equation 4.16) 
and study the pulse slippage effects. The desynchronism, as we have seen earlier in 
Chapter II, measures the shortening of the resonator cavity length by AS compared to the 
slippage distance NA. Specifically, the desynchronism between the electron pulse and 
the optical pulse is varied from d =0 to d =0.3. 

The steady-state efficiency versus the desynchronism for a step-tapered undulator 
is plotted in Figure 30. The required efficiency to achieve the 100 kW is indicated by a 


dashed line at n =0.7%. From this figure, it is obvious that we obtain the target 


efficiency and power with the conventional undulator with no taper (A=0), or witha 
negative step-taper undulator of A=-a. The desynchronism is within the interval 
d =0.03 to 0.07, with the best results for d =0.04. For the first case (A=0), the 
efficiency reaches the value of 0.8%. For the negative taper (A =~—1 ), the efficiency is 
increased to 0.86% which corresponds to an output power of 116 kW. For larger absolute 
values of tapering or for values of desynchronism out the interval 0.03 <d <0.07, the 


efficiency remains below the target value. 


The same figure verifies that at exact synchronism between the electron and the 
optical pulses (d =() there is no output power. Increasing d by a small amount, there is a 
sharp rise in the output power until the peak point of d =0.04 for these parameters. 
Increasing d beyond this point, the power slowly drops off, until eventually there is again 


no output power. 
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Figure 30. Efficiency 7 Versus Desynchronism d for Step-taper with Q, = 4.2 


Figure 31 presents a simulation output for A=-—2 with a desynchronism value of 


d =().04. The lower-left window shows the parabolic electron pulse position at tT =0 
and t=1. The middle-left window corresponds to the field ja(zn)| evolution over 
n=500 passes with the final pulse shape la(z)| shown in the upper left window. The 
color intensity scale is located on the top of this window and represents the optical field 
amplitude with the largest value of le(z)| =15.69 shown as red and zero as blue. The 


white contour line corresponds to a field amplitude of 50% of its maximum value. 
The middle lower window shows the weak field gain spectrum G(v). The 
window above that represents the evolution of the optical power spectrum P(v,n) over 


n=500 passes with the final power spectrum P(v) shown on the top window. The 


triangular tick mark indicates the initial resonance wavelength, and the rectangular tick 


mark represents the center of the final spectrum. 
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Figure 31. Simulation output for inverse step-tapered undulator. 


At the lower right window we have the dimensionless optical power evolution 
plotted after each pass over 500 passes. The optical power reaches a steady-state after 
n=50 passes. The middle-right window shows the evolution of the electron spectrum 
f (vn) over n=500 passes with the final electron spectrum f(v) shown on the 
window above it. Once more, the triangular tick mark indicates the initial resonance 


phase velocity and the rectangular tick mark represents the center of the final spectrum. 
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Figure 32. Energy Spread Ay /y Versus Desynchronism for Step-tapered Undulator. 


The slight improvement obtained with a small negative taper leads to a small 
increase of the induced energy spread Ay /y (Equation 4.14). Figure 32 plots the energy 
spread versus the desynchronism for various values of step-taper. The maximum values 
of the energy spread are found for conditions where the maximum efficiency occurs. For 
no taper (A=0) and for A=-a, the energy spread becomes Ay/y =3.5% and 
Ay /Y =3.8% respectively. Both values are well below the limit of 15% required for 


electron beam recirculation. 


We considered the possibility of decreasing undulator losses to improve FEL 


performance. The small values of the induced energy spread allow increasing the quality 
factor Q, from 4.2 to 10. The new resonator cavity has only 10% losses per pass leading 
to larger efficiency values that exceed 1% for some values of A and d. Specifically, in 
the A=—t case, the efficiency increases from 0.9% to 1.4%. With larger Q, only for 
values of A=+4m7 and for d >1.2 did the FEL failed to exceed the efficiency limit of 


0.7% [Ref. 16]. 
qZ 


In conclusion, the simulation results of the upgraded to 100 kW TJNAF FEL 
show that there is a slight improvement in efficiency with a small negative taper and the 
induced energy spread stays well below the 15% limit for the parameters investigated. In 
contrast, previously published results showed as much as a 75% efficiency enhancement 
with a negative step-tapered undulator based on simulations and experimental results 
from CLIO and FELIX [Ref. 16]. The results of our research show that the step-taper 
cannot always be as effective as found earlier, and that an FEL must be far into strong 


field saturation before tapering can extend the saturation limit. 


In the next Chapter the solid-state laser is studied. The physics behind the solid- 
state lasers is first discussed. Accordingly, their potential to be developed as directed 


energy weapon systems is presented. 
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VI. SOLID STATE LASER 


A. INTRODUCTION 

The solid-state laser is an optical oscillator. It consists of a power supply, a 
resonant optical amplifier and a feedback loop, as shown in Figure 33. An input signal is 
essential to starting the lasing process. The input is amplified and the output is fed back 
to the input with matching phase, where it undergoes further amplification. This process 
continues until the amplifier gain saturates and the system reaches a steady state with an 


output signal at the resonant frequency of the amplifier. 


Feedback 


Power supply 


Figure 33. An Optical Oscillator. (After [Ref. 17]). 





An optical amplifier is a device that increases the amplitude of an optical field. If 
the amplification takes place while maintaining the optical field’s phase, the amplifier is 
coherent and is the result of a process called stimulated emission of radiation. In this 
process, atoms, or ions at an upper energy level, are stimulated by photons to emit their 
excess energy to the input optical field coherently. Thus, the photon of the input light 
stimulates the atom in an upper energy level to undergo a transition to a lower energy 
level and emit an identical photon at the same frequency, phase, and polarization, 
propagating in the same direction. The process is continuously repeated, resulting in 
coherent light amplification. The frequency of the amplified coherent light is determined 


by the energy level difference of the amplifier material. 
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An essential ingredient for achieving laser amplification is the presence of a 
greater number of atoms in the upper energy level than in the lower level of the amplifier 
material. This is called a population inversion and is induced by the power supply, which 
excites the atoms to the higher energy level. This is a pumped active medium. The 
feedback is achieved by placing the gain medium in an optical cavity with two reflecting 
mirrors. The output is accomplished by making one of the mirrors partially transmitting 
or by constructing an unstable optical cavity. The modes of oscillation in the cavity or the 
frequency selection of the light amplification depend on the amplifier and on the 


resonator. The components of a laser are shown in Figure 34. 





ae 
Transmitted Mirror 


“ee 


Figure 34. Basic Components of a Laser. 


The solid-state lasers use solid-state material as a gain medium that can be 
electrically or optically pumped. The most common types are the neodymium doped 
yttrium-aluminum garnet laser amplifier (Nd:YAG), the ruby (Cr *:AbO3) laser 
amplifier, and the erbium-doped (Er**:Silica Fiber) laser amplifier. The lasing 
frequencies extend from the visible range to the mid-infrared while the average output 
power can be scaled to high powers (~100 kW). 

B. BASIC SOLID STATE LASER PHYSICS 

1, Radiative Processes 

Quantum theory considers the discrete energy levels in atoms and ions. The laser 
is a quantum device that is based on the interaction between the photons and the atoms 


placed in a resonator. 
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The atomic transition between two discrete energy levels, E; and E>, results in a 
radiation of a frequency f,, where hf, = E2 — E;. In a laser, there are three radiative 


processes possible: spontaneous emission, absorption, and stimulated emission. 


Spontaneous emission is the process where atoms that are initially in an upper 
energy level E2 decay spontaneously to a lower state EF), releasing their excess of energy 
in the form of photons. This transition is independent of the number of photons that 
already exist in the same mode. If the population density in state 2 is N2, then the decay 
rate of this state will be 


dN, 
dt 





=-A,N,, (6.1) 


where Az; is known as the Einstein’s A coefficient. The A coefficient represents the 
number of transitions/sec, and is equal to the inverse time constant t,, known as the 


spontaneous emission lifetime of the transition from state 2 to state 1. The inverse of the 
spontaneous emission lifetime expresses the probability density of spontaneous emission 


[Ref. 17]. 


Absorption is the process where an atom in the lower state 1 absorbs a photon 
from the incident field and transits to a higher energy state 2. The rate of this process 
depends on the absorbing material and the incident field, 


dN. dN 
- = B,N,p(f)=-—., (6.2) 
t dt 





where N77, N2 are the population densities in states 1, 2, respectively, B;2 is known as the 
Einstein’s B coefficient, and p ( ev ) is the spectral energy density of the electromagnetic 
field inside the cavity which is given by 


suf? hf _8thf>? 1 


== ES (6.3) 
p Co ole 4 = oh 
h 
<E>=—5 if : (6.4) 
eke | 
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where k is the Boltzmann’s constant, h is the Planck’s constant, c is the light’s 


propagation speed, and <E> is the average energy per mode in thermal equilibrium. 


Stimulated emission is a reverse process of absorption, where the atom is in the 
upper energy level and is stimulated to emit a photon into the same mode of the incident 
field, resulting in coherent photon amplification. The rate of the stimulated emission 


process is given by 








dN, dN, 
=—B,,N )=-—, 6.5) 
dt 21°" 2 p (f dt ( 
These three radiative processes are shown in Figure 35. 
N2 
E, hf 
hf hf 4\S\AW2hf 
A B B 
21 L\S\V N/a 12 7\S\V 21 AW 
(a) (b) (c) 
Figure 35. (a) Spontaneous Emission (b) Absorption (c) Stimulated Emission. 


The total rate of change of the population in state 2 or in state 1 of the material as a result 
of the radiative process is the sum of Equations (6.1), (6.2), and (6.5), 
dN, dN, 


=—A,,N,+B,,N.p(f)—-B,,N,p(f) =- : (6.6) 
dt dt 








At thermal equilibrium, the total rate must be zero, and Equation (6.6) provides the ratio 
of the population of the two states, 


N, A,t+ Bef) 
The ratio of the population of the two states according to the Boltzmann probability 


distribution is given by, 


hf 
Ma 2 8a, (6.8) 
Nie Gi 
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The degeneracy parameters g,,g, represent the number of states corresponding to 
the energy levels E; and E> respectively. If Equations (6.7) and (6.8) are combined, 
solved for (jf) and set that value equal to Equation (6.3), the relation between the 


Einstein’s coefficients will be given as, 





By _ Ba (6.9a) 
By, 8 

Sr hf 
ee (6.9b) 
B,, Cc 


2. Spectral Line Widths 

In reality, a system consists of many atoms and there is an inherent uncertainty in 
state energies due to finite upper state lifetimes, Doppler shift of moving atoms, and 
collisions between the atoms. This means that there is a distribution of energies around 
the mean values of state energies, as shown in Figure 36a, and the atomic emission or 
absorption has a relative narrow spectral distribution. Thus, the system emits or absorbs 


radiation in a narrow band of frequencies of width Af . This relative distribution is called 
the lineshape function g(f) and represents the relative strength of the atomic interaction 


with photons of different frequencies. The lineshape function is centered at central 


frequency fo and is normalized so that 
| s( fap '=1. (6.10) 
0 


Therefore, transitions are most likely to occur for photons of frequencies f = f,. 
The width of the function ral I ) is the transition linewidth Af and is defined as the 


width of g ( a ) at its half maximum value (FWHM). The function g ( - ) is shown in 


Figure 36b. 
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Figure 36. (a) Energy Level Diagram ___(b) Lineshape Function g/(f). 


The atoms in the material are not dentical, and they have different lineshape 
functions and even different central frequencies. In this case, the composite lineshape is 
defined, which is the average lineshape with its peak shifted from either central 


frequency. 
Taking into account finite Inewidths, the radiative rate Equation (6.6) for N2 
becomes 


dN. T T T 
== AN, | Pdf + BN, | PCE 8 df BN. | PCP 8S df". (6-11) 





Assume that the spectral width p(f) is very small compared to g(f).So p(f) 


can be approximately represented as a 6 —function, p(f)=p,5(f'— f). Then the rate 


equation becomes 








dN 
7 = —A,,N,+B,N,P ;-8(f)— By Np -a(f). (6.12) 
Combining the last equation with Equations (6.9a) and (6.9b), we get 
dN. I, 
2 Aa -| Ne Say heen, (6.13) 
dt & Jaf 


where I, =P; (c/n) is the light’s intensity, n is the index of refraction of the material, 


and © ( if ) is the stimulated emission cross section area, 
x? 
O(f)=Ay——> af). (6.14) 
Sn 
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3. Laser Amplification 

Light transmitted through a material in thermal equilibrium is attenuated due to 
absorption by the large population in the lower energy level. With the population 
inversion where Nz > N;, a nonequilibrium situation is achieved and therefore laser 


amplification is possible due to stimulated and spontaneous emission. 
Let J, be the incident light intensity to the amplifier material of length dz. The 


incremental emitted intensity d/, from within that length will be 


Snip S, (6.15) 

The radiation due to spontaneous emission is spread uniformly into a solid angle 
of 42 steradians. Therefore, of the spontaneously emitted photons, only those in a solid 
angle dQ are in the correct direction and only half of that fraction has the proper 


polarization. Thus, Equation (6.15) can be rewritten, 


f 2 j 
—=hf | —-——N,,A,,+| N, -—N, /-o : 6.16 
dz i} An 24721 f 2 ; i a) ( ) 
If considering only the coherent amplification and ignore the spontaneous emission, we 
have 
dl, 82 
=~ =|N,-=N, OCP, SYP,» (6.17) 
dz 81 


where y(f ) is the gain coefficient and represents the net gain in the light intensity per 


unit length of the medium. The solution of Equation (6.17) is the exponentially increasing 


function, 


I(z)=1, Oe. (6.18) 


It is obvious that the gain coefficient is y(f) proportional to the population 


difference of the two energy levels. With population inversion, the gain coefficient is 


positive and there is light amplification. In the absence of population inversion, the gain 
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coefficient is negative and the medium will attenuate rather than amplify the light 


traveling though it. 


For a total length L of the material, the overall gain G(f ) of the laser amplifier is 

the ratio of the light intensity at the output to the light intensity at the input, so that, 
G(f)=e"™, (6.19) 

From Equations (6.14), (6.17), and (6.19), it is obvious that the gain coefficient 
y(f) and the gain G(f) are proportional to the lineshape function g(/) and are 
dependent on the frequency of the incident light f. Therefore, the laser amplifier is a 
resonant device, with a resonant frequency and bandwidth determined by the lineshape 
function g ( 7c ) 

4. Gain Saturation 

Consider the energy diagram of a two level system with a ground state energy 
level Eo as shown in Figure 37. The levels 1 and 2 have overall lifetimes Tt, and T, 
respectively, for transitions to lower levels. If R; and Rp are the pumping rates of the 


levels 1 and 2 respectively, then the basic Einstein’s rate Equation (6.13) will be for level 


2, 








dN. 
: =R,(@) -A,,N,-W(f)[N,— N|-AnN. (6.20) 
and for level 1, 
dN, 
=R(t)+A,N,+W(f)[N-N,J-AN,, (6.21) 


where now g,=g, and w( a: ) =O ( he ) I, / hf . Since the rates of transition are inversely 


proportional to the decay times and A209 is usually small, the last two equations can be 


rewritten as, 





= = R)-=2-W(P[N Ni) (6.22) 
Ai = RB ewip[y,-W |, (6.23) 
dt T, Tv, 


where T5'=135) +75. 
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Eo 
Figure 37. Two Energy Level 1 and 2 System with Ground State Eo. 


For steady-state conditions of the populations of the levels 1 and 2, the last two 
equations for the population difference AN can be solved 


Rt,-(K + RM, 


AN = N,- N= 
14+7,W(f) 


(6.24) 


At t=0 there is no field, W(f) =0 and the population difference at this time has 


its maximum value AN(0)=R,t,—-(R,+R,)t,, which represent the steady-state 


population difference in the absence of amplifier radiation. Therefore, in the presence of 


amplifier radiation the steady state population will become 


AN (0) 


I,” 
144 
I 


Ss 





N,-N,= (6.25) 


where J/, =hf i o(f)t, is the saturation intensity. When / y =/,, the population 
difference is one half the maximum value AN (0) . As previously seen from the Equation 


(6.17), the gain coefficient is proportional to the population difference. Thus, the 


saturated gain coefficient (f ) can be defined as 


v(f) 1D (6.26) 


14+—4 
I 


Ss 


where y,(f)=AN(O)o(f). It is obvious that the gain coefficient is a decreasing 


function of the intensity 7,. The saturation intensity is defined as the intensity that 


i 
reduces the gain coefficient by a factor of 2. 
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3 Three and Four Level Lasers 

A population inversion arises when more atoms are in a higher state of excitation 
than in a lower state. This can occur when the upper state has a longer lifetime. A large 
number of atoms can be excited and held in this state leaving an almost empty state 
below it. The atoms can stay in this metastable state while the population inversion is 


built up. 


In reality, it is not possible to have lasing based on absorption and emission 
between only two energy states. For any pair of states, the rate at which the upper level is 
populated by absorption equals the rate at which atoms decay to the lower level by 
stimulated emission. Three or four energy levels schemes are used in practice to achieve 


a population inversion. 


A three level system is shown in Figure 38. The atoms are pumped to the upper, 
short-lived level £2), which is called the pump level. A rapid spontaneous decay occurs 
from the pump level to the metastable level, which lies between the pump level and the 


ground state. Laser emission occurs between the metastable level and the ground state. 


An improvement on this scheme is obtained with a four level structure which is 
shown in Figure 39. Laser transition takes place between the third and the second excited 
states. The depopulation of the highest and the lower laser level is rapid in order to ensure 


that always a population inversion exists. 
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Figure 38. Three Energy Level System. 
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Figure 39. Four Energy Level System. 


6. Optical Resonant Cavity 

It has been demonstrated that to obtain light amplification from a given medium it 
is necessary to establish a population inversion between a chosen pair of energy levels 
and promote stimulated emission. This by itself is not enough to sustain laser action. The 
lasing medium must be bounded between two mirrors, of which, one is totally reflecting 
and the other is partially reflecting. This forms an optical resonant cavity where the light 
oscillates between the mirrors. The cavity provides a method of optical feedback where 
the light is reflected several times between the mirrors to enable further stimulated 


emission. 


A cavity of length S$ supports only the modes such that S corresponds to integral 
multiples of half a wavelength. The resonator mode spacing due to the round trip phase 
shift is c/2S, where c is the speed of light in the medium. Furthermore, laser amplification 
can occur only for the modes of frequencies for which the signal gain coefficient y, (f) 
is greater than the loss coefficient a, in the cavity. The coefficient loss expresses the total 


number of photons lost per unit length, and can be defined as [Ref. 17] 


a,=—., (6.27) 


where T,, represents the photon lifetime. After a short time of laser operation, the gain 


saturates to y(f) and only a single or a few modes close to the central mode survive. 


Finally, all the emitted photons have the same wavelength, phase, and amplitude. The 
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process of mode competition is shown in Figure 40. The plot on the left shows that, when 
the hser is turned on, all the resonator modal frequencies for which the gain coefficient 
is greater than the loss coefficient (presented as blue lines) begin to grow, with the central 
mode growing at the highest rate. The modal frequencies for which the loss exceeds the 
gain are indicated as red lines. When the gain saturates, only a single mode survives as 


shown on the right plot. 


The laser can support either a continuous wave (CW) beam or a pulsed beam. The 
pulsed laser beam is achieved either by controlling the pumping rate (gain switching) or 
by controlling the cavity losses (Q-switching). In a gain switching laser, the pumping rate 


is turned on and off periodically. In a Q-switching laser, the loss coefficient a, is 


switched between a large value during the off time and a small value during the on time. 
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Cavity 1 \ 
Losses 
, — coy “a heerak:, inde le, el 
/ \VWf) 
/ \ 


nan | 
Allowed f Mode of f 
modes oscillation 
Figure 40. Mode Competition. 


C. SOLID STATE LASER MATERIALS 

This section will discuss some of the large number of solid sate materials that 
have the potential to be scaled to high average output power for military applications. 

1, Nd: Lasers 

Neodymium-doped laser materials are the most commonly used as they have a 
combination of properties uniquely favorable for laser operation. They are excited near 


808 nm and emit near 1.06 um. 
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The neodymiumdoped yttrium aluminum garnet (Nd:YAG) laser, because of its 
high gain, is the most important. The YAG material is hard, of good optical quality, and 
has high thermal conductivity. 


The Nd:YAG laser is a four level system. Under normal operating conditions at 
room temperature, the energy levels allow a transition at 1.0641 um _ wavelength. 
However, it is possible to obtain lasing at other wavelengths either by inserting etalons or 
prisms in the resonator, or by designing a resonant reflector as an output mirror. A design 
and operation of a high power Nd:YAG laser oscillating at 1.3 um has been reported. In 
the long pulse mode the laser has achieved an output up to 8 J per pulse at 5 Hz and up to 


165 W of laser output at 50 Hz [Ref. 18]. 


Lithium yttrium fluoride lasing material has been improved and is suitable for 
many applications. The Nd: YLF has a relatively large thermal conductivity which allows 
efficient heat extraction, but less than Nd:YAG. Its excited state lifetime is twice as long 
as Nd:YAG. Since lifetime is a measure of the energy storage capability in Q-switched 
operation, Nd: YLF is currently being used in Q-switched laser amplifiers. The stimulated 
emission cross section area (6.14) for YLF is also half that of YAG. Depending on the 


polarization of the light, two wavelengths are obtained, 1.05 um and 1.3 um. 


Another material is the neodymium-doped yttrium vanadate (Nd:YVO4), which 
has a large stimulated emission cross section area, five times higher than Nd:YAG, and a 
strong broadband absorption. The achieved high gain combined with the strong 
absorption allows us to use a crystal only a few mm in length. Pump absorption is 
polarization dependent. The strongest absorption occurs for pump light polarized in the 
same direction as the light propagation. However, the thermal conductivity is only half as 
high as Nd:YAG and somewhat lower than Nd:YLF. The excited state lifetime is also 
shorter than Nd:YAG. Therefore, the large cross section is partially offset by its shorter 
lifetime. 

De Yb: Lasers 

Ytterbium-doped materials (Yb: YAG, Yb:FAP) have been known for years but 


were never considered interesting due to their low absorption and emission cross- 
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sections. Moreover, they appear to have a quastthree-energy-level nature, since the 
terminal state can be thermally populated. By comparison, the separation between the 
terminal laser level and the ground state in Nd:YAG at room temperature is ten times 
bigger than in Yb:YAG. The crystal of Yb: YAG has only a single absorption line at 942 
nm, which requires extremely high pump irradiance with conventional pump sources in 
the visible band to achieve lasing threshold. However, since bright diode-pump arrays of 


InGaAs are available, which emit at 942 nm, the Yb: YAG has become attractive. 


The Yb:YAG is a quast three-energy- level laser amplifier and the wavelength of 
oscillation is 1.03um. Its performance is strongly dependent on the operating 
temperature and can be improved by cooling the crystal. In this way, the thermal 
population is reduced and the stimulated emission is increased. To accomplish this, the 
crystal is mounted on a heat sink, which could be cooled to low temperatures. Thus, laser 
efficiency and the cost of the pump source can be traded for overall system efficiency 


that includes power and refrigeration requirements. 


Despite these difficulties, Yb: YAG has attractive properties and the potential for 
reaching higher output power than other materials. The thermal load generated in a laser 
medium is a great concern for high power scaling. A reduction in heat load can lead to 
higher power systems with good beam quality. A unique advantage of Yb:YAG is that, 
when pumped with an InGaAs pump source, it produces a smaller amount of crystal 
heating than any other known material. The required pump radiation in this material is 
only about one-third compared to Nd:YAG. Another unique advantage is that it has a 
broad absorption line at 942 nm, which is ten times broader than the 809 nm absorption 
line in Nd:YAG. It has also a long lifetime which reduces the number of pump diodes 
required for a given energy per pulse output. As will be seen, the pump diodes are power 
limited, and the long lifetime permits a long pump pulse, which generates high power 
output. In addition, the InGaAs family of diodes offers greater lifetime and higher 


reliability than any other diode family. 


The obvious advantages of very low thermal load, broad absorption at the InGaAs 


wavelength, and long lifetime of Yb:YAG have raised the prospect of reaching high 
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average power. This feature has motivated the research and development of ytterbium- 
doped materials. Figure 41 indicates the progress in development of the Yb:YAG laser 


over the past few years. 
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D. SOLID STATE LASER PUMP SOURCES 

The pump source of a solid-state laser emits radiation in a spectral region close to 
the absorption bands of the lasing material. The most efficient laser pump radiates at 
wavelengths that match the useful absorption band of the lasing medium. The pump 
source is supplied by electrical energy, either continuous or pulsed, which the pump 
converts into optical radiation. An optical system may be used to couple the radiation 


from the pump source to the laser medium. 


Today, the pump sources that are of practical interest are flashlamps, cw arc 
lamps, and laser diodes. Coherent pump sources, such as laser diodes, permit 
concentration of the pump radiation into the active material with few losses. For lamp 
sources, only a fraction of the omnidirectional radiation can be concentrated into the 
active medium. Typical overall electrical efficiency for a lamp pumped laser is 1%, for an 


arc lamp about 3%, and for a diode pumped lasers 10-15% [Ref. 18]. Until recently, the 
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low power output, low packing density, and extremely high cost of diode lasers prevented 
any potential application of them. The rapid advances in laser diode technology have 


made the diode pump sources the most efficient method of pump. 


Actually, flashlamps have higher radiation output to electrical input efficiency 
(70%) compared to laser diodes (30%). However, the diode pump is the most efficient 
way to produce population inversion in the laser material because there is an excellent 
spectral match between the output from the diode pump and the laser medium. This 
reduces the amount of heat deposited in the laser material and therefore leads to better 
beam quality. Since the output of a diode array is directional, the radiation can be 
transferred to the laser medium with little loss. Moreover, the diodes are coherent devices 
and their beams can be focused and adjusted to match the resonator modes of the laser, 


which in turn, leads to a high brightness laser output. 


The diode pump can be configured in two ways: the side pump and the end pump. 
In the side pump geometry, the diode arrays are placed along the laser material while in 
the end pump geometry, the pump radiation is focused into the laser material collinear 


with the resonator mode. 


Laser diode arrays have lifetimes on the order of 10* hours in cw operation or 10? 
shots in the pulsed mode [Ref. 2]. This is a critical factor as it influences the system 
lifetime and reliability. Moreover, new laser diodes, such as the InGaAs family, enable 
the use of laser materials such as Yb:YAG, which seems to be very promising for 


reaching high average power output. 


Diode laser arrays have a major impact on the entire laser system’s performance, 
cost, and reliability. Primary diode array characteristics are output power and duty factor, 
operating temperature and cooling, spatial and spectral characteristics, reliability, 
lifetime, and cost. Diode array performance is affected by the characteristics of the diode 
bars themselves, by the coolers on which they are mounted, and by the laser environment. 
Technological advances in these areas are required in order for solid-state lasers to realize 


their potential for directed energy weapons applications. 
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E. THERMO-OPTIC EFFECTS AND HEAT REMOVAL 

The pumping process of the gain medium in the solid-state laser is associated with 
the generation of heat, causing various thermal effects. The thermal effects are considered 
the primary design problem for high average output power systems. The generation of 
heat in the laser medium during the pumping process is caused by: 


e Quantum defect heating, which is the energy lost as heat to the lattice 
since there is an energy difference between the pump band and the upper 
laser level. In a typical laser system, the transition frequencies between 
these two levels fall within the frequency range of the vibration spectrum 
of the lattice. Therefore, these transitions are accomplished by direct 
nonradiative decay (emitting phonons), which leads to the generation of 
heat. 


e Similarly, the energy difference between the lower energy level and the 
ground state. 


e The fact that the quantum efficiency of the emission processes in the laser 
is less than unity. The quantum efficiency relates the number of 
independent ionic transitions to the number of photons generated by this 
process, and depends solely on the position of the energy levels [Ref. 18]. 


e The bw pump efficiency especially when arc lamps or flashlamps are 
used. 


The heat from the pump radiation and the surface cooling of the laser medium 
cause a nonuniform temperature distribution across the laser medium. This temperature 
variation leads to a variation of the index of refraction and thermal effects such as 
thermal lensing and stress birefringence (induced birefringence that leads to beam 
depolarization). Both of these affect the laser beam quality. Moreover, thermal load 
causes stress on the material. If the temperature of the material exceeds its tensile 


strength, stress fractures occur in the medium. 


A method to correct for thermal effects on the laser beam is the phase conjugation 
technique [Ref. 18]. Suppose a uniform flat wavefront which is distorted in the amplifier 
medium. An ordinary cavity mirror inverts the distortion as it reflects the beam, keeping 
the distortion fixed with respect to the propagation direction. With a second pass through 
the amplifier, the distortion is doubled. Using a phase-conjugate mirror, the distortion is 
reversed relative to the wave propagation direction. Thus, the same region of the 
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amplifier that created the distortion compensates it during the second pass through the 


amplifier. 


Phase conjugation has been used successfully to correct the aberrations induced to 
the laser beam by thermal and stress effects for multi-pass, single solid-state amplifier 
configurations. In this area, TRW has already demonstrated clean-up of a 200 W average 
power beam [Ref. 2]. Still, there is progress to be made in using phase conjugation to 
combine multiple solid state amplifiers in phase in order to achieve high average power 
sources. Until now, Lawrence Livermore National Laboratory (LLNL) and Raytheon 
have independently demonstrated phase combining of four or more amplifier chains [Ref. 


Z|: 


There are several methods to remove the unwanted heat from the laser system 
such as liquid cooling, air or gas cooling, and conductive cooling. The cooling techniques 
remove the heat generated in the laser medium, pump source and laser cavity. In liquid 
cooling, the coolant flow is under pressure over the lasing medium and the pump source 
surface. Air or gas cooling is usually used for low average power lasers. In conductive 
cooling, a heat sink is mounted directly to the laser material and the heat is removed 


through conductivity. 


Cooling is a very important process for laser operation, especially for a MW level 
solid-state laser. The role of the cooling process is obvious when considering the energy 
flow inside the system for a single engagement. An engagement using this type of 
weapon system requires an output of 1 megajoule (MJ) energy from the beam director. A 
typical laser medium is only 40% efficient, which means that the required energy from 
the pump source is 2.5 MJ. The pump source itself, as mentioned in the previous section, 
is expected to be 50% efficient. Thus, the electrical power supply should provide a total 
energy of 5 MJ. An energy flow diagram from the power supply to the laser output is 
shown in the Figure 42. From that figure, it can be seen that the cooling process should 


absorb 4 MJ unwanted heat from the weapon system components. 
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Figure 42. Energy Flow for a Single Engagement.(After [Ref. 2]). 
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F. HIGH ENERGY SOLID STATE LASER DESIGNS 

A considerable amount of progress has been made over the last 15 years in scaling 
solid-state lasers to higher powers. Kilowatt average power levels are now routine and 
10 kW output will soon be demonstrated. Systems at the 100 kW level are being 
designed today, and conceptual designs for MW level devices have been developed by 


multiple industrial sources. 


Solid-state lasers with an average power of 0.1-10 kW are currently in use for 
many military applications. In addition to their typical role as designators and 
illuminators, lasers of this category can be used for countermeasures and sensor denial, 


but they are not generally considered directed energy weapon systems. 


Solid-state lasers with an average power of 10-100 kW can place a sufficient 
amount of energy density on military targets to cause their destruction. However, at this 
power level, the lethal engagement of the target is limited to close ranges, especially 
when the laser weapon system has to accomplish its mission in the maritime 


environment. 


Lasers with an average power of 100-500 kW are typically associated with high- 
energy laser weapons. The Theater High Energy Laser (THEL), developed for defense 
against Katyusha rockets, the Heat Capacity Laser being developed for use in the Army’s 
Enhanced Area Air Defense mission, and the proposed Advanced Tactical Laser (ATL) 


are in this class. These types of lasers produce lethal energy fluence on a target during 
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short engagement times, at ranges from 3 to 50 km depending upon the target, 


engagement altitude and local atmospheric conditions. 


Lasers with an average power of over 500 kW are of interest to the Navy for point 
defense against missiles, and for strategic applications for ballistic missile defense over 


long distances. The Airborne Laser and the Space-Based Laser fall into this class. 


High power solid-state lasers can be grouped into four classes: diode lasers, fiber 
lasers, continuously cooled lasers and heat capacity lasers. 

1. Diode Lasers 

Diode lasers use a variety of semiconductor compounds and can be either 
electrically or optically pumped. Lasing has been demonstrated throughout the visible 
range and extends well into the mid-infrared wavelengths, with the highest achieved 


average powers in the near infrared band. 


Single diode lasers, using special semiconductor architectures, have achieved 
power levels of 1 W in CW mode. In order to achieve high average power output, 
multiple diode amplifier chains must be used. Diode arrays have demonstrated more than 


1 kW/cn?’ output intensity. 


Incoherent beam combining (IBC) has been used successfully to couple multiple 
diode lasers. An experiment has demonstrated an output of 5 W with acceptable beam 
quality. Current research at the Air Force Research Laboratory (AFRL) is trying to scale 
the power output to multiple tens of Watts with good beam quality, and a laboratory 
demonstration is expected within the next years. However, this technology is still 
immature, and it remains to be seen whether it is a scaleable technology. The potential 
exists for achieving 60% efficiency with excellent beam quality by combining multiple 


diode emitters. 


The Phase Integrated Laser Optics Technology (PILOT) program, conducted at 
AFRL during the late 1980s and early 1990s, demonstrated coherent phasing of 900 
diode elements, but the total output power was less than 10 W. The PILOT research and 
development effort was discontinued, but similar phasing techniques may be applied to 


current fiber laser designs [Ref. 2]. 
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2. Fiber Lasers 

Low power versions of these lasers are commonly used in fiber optic 
communication systems. High power fiber lasers have become enabled by the 
development of fibers that have a doped core separate from the cladding. This 
configuration provides a much larger aperture coupling in the pump light, so higher 
pumping levels are achievable, which leads to higher average power output. Continuous 
wave cladding-pumped fiber lasers have achieved an average power exceeding 120 W 


and a conversion efficiency of diode pump light to output laser radiation exceeding 60%. 


Figure 43 displays the different configurations of the fiber lasers. On the left, a 
conventional fiber amplifier is shown where an end diode pump process is configured. 
Separating the core from the cladding, a side pump process is accomplished as shown on 
the right. The fiber core is indicated as a red cylinder while the cladding is indicated as a 


yellow cylinder. 





Diode-Pump 
Light 
Core 
Polymer Cladding Cladding 
Conventional Fiber Laser Cladding-Pumped Fiber Laser 
Figure 43. Conventional and Cladding- Pumped Fiber Laser Configurations.(From 
[Ref. 2]). 


The potential of fiber lasers for directed energy weapons applications lies in two 
unique features: 

e The central core can be designed as a single-transverse- mode waveguide 

e the pump aperture of any individual fiber is on the order of 200 microns. 

The single-transverse mode operation is achieved when the central core is limited 
to less than 10 microns. These small core fibers will produce diffraction limited output, 
but are limited in power to approximately 100 W by nonlinear phenomena in the fiber 


such as Stimulated Raman and Brilluouin scattering. New methods of photonically 
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engineered fibers with a radially varying refractive index profile within the core offer the 


possibility of a larger single-mode. 


By combining one thousand fibers, each carrying an average power of 100 W, a 
100 kW laser system with a cross sectional aperture less than 1 cn? in area could be 
produced. However, at this point, it is necessary to deal again with the phase combining 
problem. Each individual fiber is an independent laser and, as a result, the output of the 
individual fibers adds incoherently. Thus, the fiber array indeed produces 100 kW output 
power, but in a highly divergent beam, with a time-dependent, random coherence. The 
brightness of the light beam of incoherently combined arrays will be limited compared to 
coherently combined arrays, but may still be sufficiently large to be valuable for some 


applications. 


In Figure 44, an array of independent fibers produces a highly divergent beam, 
which is also a characteristic of a single fiber output with a very small aperture and a 
time-dependent speckle pattern. In the same figure, the desired output beam of the 
coherently phase combined fibers can be seen. In this case, the beam is free of speckle 


and exhibits the divergence characteristics of the effective aperture of the entire array. 


The promise of fiber lasers for high-energy applications can be realized if the 
technology to coherently phase combine the output of individual fibers can be developed. 
This is an active research area in many laboratories. The state-of-the-art achievement is 
two low power fibers successfully phase-combined together producing a coherent output 
of 100 mW [Ref. 2]. However, in order to accomplish the goal of a MW output, phase 


combining of a large number of fiber lasers is required. 





Independent Fiber Lasers Coherently phased fibers 


Figure 44. A Beam from an Array of Independent Fibers and a Beam from 
Coherently Phased Fibers.(From [Ref. 2]). 
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3. Continuously-Cooled Lasers 

The continuously cooled lasers have already achieved over 1 kW output power. 
During their operation, the cooling process of the laser medium occurs simultaneously. 
Both the amplifier medium and the cooling method were discussed in previous sections. 
The most promising configuration, as demonstrated, is that using ytterbium-doped 
amplifier materials pumped by the InGaAs diodes. The phase combining of a large 
number of these apertures has the potential to achieve the desired output power for solid- 
state laser weapon. 

4. Heat Capacity Lasers 

The heat capacity laser has already achieved an output of 1.4 kW power at LLNL, 
and a power scaling up to 100 kW is expected in the near future. In contrast to the 
continuously cooled lasers, lasing and cooling occur separately in the heat capacity laser. 
With this method, there are two cycles of the laser operation: the lasing cycle followed by 
the cooling cycle. After each laser emission, which means a target engagement for the 
weapon design system, the cooling cycle removes the heat that has been deposited in the 
system. Typically, with current cooling technology, after an engagement of a couple of 


seconds, a cooling period of | to 5 minutes follows. 


The heat capacity laser has been chosen for the Army’s tactical high-energy laser 
(HELSTAR). The program is based on the solid-state heat capacity laser developed at 
LLNL. This type of laser operates in pulsed mode with a typical width of several hundred 
microseconds and a repetition rate of several hundred Hertz. In order to achieve the goal 
of 100 kW output power, a diode pumped design will replace the current flashlamp 
pumped laser system and Nd:GGG crystalline will replace the current amplifier material 
of Nd:glass. 

G. CONCLUSIONS 

In this Chapter the solid-state laser was studied. This laser also offers significant 
potential for a directed energy weapon system since it has already demonstrated a kW 
level output power. In order for solid-state lasers to realize their potential for directed 


energy application, promising laser designs, such as the fiber, continuously cooled and 


oF 


heat capacity lasers must develop independent technologies dealing with the laser 


materials, pump process, thermal management and beam control. 


Finally, the following Chapter summarizes the results presented in this thesis. 
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Vil. CONCLUSIONS 


High-energy laser (HEL) systems provide military units the opportunity to 
improve their capabilities and tactics in a variety of mission areas. A naval HEL 
application is anttship missile defense at a range much greater than that offered by 
currently used weapon systems. The Navy has supported research that is expected to arm 


future ships with free electron or solid-state lasers. 


At the Naval Postgraduate School, a high power (MW level) FEL is being 
designed in collaboration with Jefferson Laboratory and Advanced Energy Systems. The 
designed FEL should have a small footprint and compact configuration in order to be 
placed onboard a Navy ship. To satisfy the size requirements, the size of the optical 
cavity is limited and the mirrors are vulnerable to the powerful optical beam. The use of a 
short Rayleigh length undulator may be the solution to the problem since it increases the 
spot size at the mirrors and leaves a small optical mode waist in the center of the cavity. 
For a small optical mode waist, a slight misalignment of the electron beam, either an 
offset or a tilt, reduces the overlap between the electrons and the optical mode in the 


interaction region, resulting in less gain and efficiency. 


In this thesis, a three-dimensional simulation is used to demonstrate that the 
electron beam shift can change the optical wavelength, and reduce the peak efficiency as 
the electron beam is offset further from the undulator axis. It is found that the required 
efficiency for the 1 MW goal (0.7%) is achieved for electron beam shifts less than 0.3 
mm, which is well beyond the experimental design tolerance of 0.01 mm. Electron beam 
focusing was also studied by varying the electron waist radius while keeping the total 
current fixed. This focusing enhances the overlap between the electron beam and the 
optical mode. It is found that the focused electron beam increases the extraction 
efficiency far beyond the desired value of 0.7%, and indicates that it may be possible to 
obtain MW operation at a lower average current. The induced electron energy spread for 


the focused beam is increased, but remains within the design goal of 15%. 
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The proposed FEL design operates in a pulse mode, producing optical 
micropulses 1 ps long at a pulse repetition frequency of 02=750 Mhz. The pulse mode 
of operation leads to frequency spread depending on the FEL pulse shape in the time 
domain. Propagation of a broad-spectrum laser pulse through the atmosphere increases 
the power in frequency components where the absorption coefficient can be high. This 
increases the total absorbed power, which may cause thermal blooming. It is found that 
thermal blooming increase is not an issue for a moderately clear atmosphere when the 


stagnation zone size remains less than 10 m. 


The TJNAF FEL can be upgraded to operate at 100 kW in the near future. It was 
shown that there is a slight improvement in efficiency with the use of a small negative 


step-tapered undulator, while the induced energy spread stays well below the 15% limit. 


In the last chapter of this thesis, the solid-state laser was studied. This laser also 
offers significant advantages in the design of the directed energy weapon system. 
However, in order for solid state lasers to realize their potential for a directed energy 
weapons application, they must be developed as high average power, high energy pulse, 
lightweight, compact systems. In this area, technological advances are required in both 
the supporting technologies and the laser itself. The heat capacity laser is about to reach 
100 kW showing tremendous progress, but this is still a laboratory device and not a field 
system. Moreover, the 100 kW average power is still not enough for the Navy’s 


requirements for a megawatt class laser design. 
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